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ABSTRACT 

The  primary  goal  of  this  thesis  was  the  development  of  an 
experimental  technique,  and  supporting  software,  for  the  ac- 
quisition and  analysis  of  unsteady  velocity  data  generated  by 
an  oscillating  airfoil.  This  research  was  in  support  of  a 
major  investigation  of  the  compressibility  effects  on  dynamic 
stall  . 

The  experimental  procedure  involved  schlieren  flow  visu- 
alization for  comparison  of  steady  and  unsteady  flow  fields, 
and  for  determination  of  parameters  for  further  study.  Laser 
Doppler  velocimetry  was  employed  for  obtaining  velocity  data 
in  the  airfoil  wake.  For  unsteady  data,  the  airfoil  was 
oscillated  in  pitch  about  its  quarter  chord. 

The  data  analysis  produced  wake  profile  plots  represent- 
ing the  flow  field  disturbed  by  the  airfoil.  Results  were 
obtained  for  steady  and  unsteady  conditions. 
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I  .   INTRODUCTION 


A.   DYNAMIC  STALL 

The  possible  benefits  and  consequences  inherent  in 
dynamic  stall  have  stimulated  experimental  investigations 
leading  to  significant  progress  in  understanding  this  complex 
phenomenon.  Much  of  the  earlier  work  in  this  area  was  car- 
ried out  by  McCroskey  using  a  two-dimensional  airfoil 
section  oscillating  in  pitch  about  its  quarter-chord.  The 
results  of  this  effort  succeeded  in  clarifying  the  features 
of  dynamic  stall  up  to  a  Mach  number  of  0.3  CRef.  11. 

Several  significant  distinctions  between  static  and 
dynamic  stall  should  be  observed.  Static  stall  is  character- 
ized by  massive  flow  separation  and  loss  of  lift  coinciding 
with  some  critical  angle  of  incidence.  The  aerodynamic 
forces  on  the  airfoil  vary  uniquely,  or  nearly  so,  with  angle 
of  attack.  Dynouriic  stall,  occuring  during  rapid  increases 
in  incidence,  results  in  a  delay  of  the  associtated  loss  of 
lift  until  well  beyond  the  static  stall  angle.  Dynamic  stall 
also  tends  to  be  considerably  more  persistant.  The  flow 
field  does  not  immediately  adjust  to  reductions  in  incidence 
angle  as  with  static  stall.  CRef.  23 

Dynamic  stall  is  characterized  by  the  shedding  and  pas- 
sage of  a  vortex  over  the  upper  surface  of  an  airfoil.  The 
production  of  this  vortex  begins  shortly  after  the   incidence 


has  exceeded  that  of  static  stall.  A  thin  layer  of  reversed 
flow  in  the  boundary  layer,  evolving  into  a  vortex  near  the 
leading  edge,  represents  the  initial  phase  of  a  dynaniic  stall 
cycle.  As  this  disturbance  moves  across  the  upper  surface  cf 
the  airfoil,  large  negative  pitching  moments,  significant 
drag  rise,  and  a  monotonic  increase  in  lift  result.  The 
aerodynamic  forces,  reaching  a  maximum  as  the  vortex  ap- 
proaches the  trailing  edge,  decrease  dramatically  with  the 
shedding  of  this  vortex.  Reattachment  of  the  flow  will  then 
typically  occur  at  some  angle  less  than  that  corresponding  to 
static  stall.  CRef.  23 

B.   THESIS  GOALS 

The  research  project  generating  this  thesis  represents  a 
small  part  of  a  major  investigation  aimed  at  developing  a 
thorough  understanding  of  the  effects  of  compressibility  and 
pitch  oscillation  on  dynamic  stall.  The  primary  goal  of  the 
research  associated  with  this  thesis  was  the  development  of 
an  experimental  procedure,  and  supporting  software,  for  the 
acquisition  and  analysis  of  unsteady  flows  with  a  laser  dop- 
pler  velocimetry  system.  The  procedure  and  software  devel- 
oped was  to  provide  an  initial  analysis  of  possible  areas  of 
interest  and  a  base  upon  which  to  build  as  the  sophistication 
of  the  experimental  methods  increased.  A  successful  attain- 
ment of  this  goal  required  obtaining  velocity  data  in  the 
wake  of  an  oscillating  airfoil  and  wake  profile  plots  pro- 
duced from  this  data. 
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I  I  .   EXPERIMENTAL  CONSIDERATIONS 

A.   EXPERIMENTAL  METHODS 

1  .   Laser  Doppler  Veloc  imetrv 

Since  its  introduction  in  1964,  laser  Doppler  vel- 
ocimetry  has  proven  to  be  a  highly  reliable  means  for  obtain- 
ing measurements  in  complex  turbulent  flows.  The  method  is 
capable  of  performing  instantaneous,  non- i ntrusi ve  velocity 
measurements  under  a  wide  range  of  conditions.  Being  rela- 
tively independent  of  fluid  properties,  accurate  measurements 
are  possible  without  the  requirement  of  calibration.  CRef.  31 

The  basic  principle  behind  laser  Doppler  velocimetry 
is  the  measurement  of  the  velocity  of  particles  travelling 
with  the  fluid.  These  particles,  when  illuminated  with  a 
focused  laser  beam,  become  sources  of  scattered  light.  The 
Doppler  shift  of  the  scattered  light  is  then  used  to  deter- 
mine the  velocity  of  the  particles  and,  depending  upon  how 
well  the  particles  are  following  the  flow,  the  fluid  velocity 
as  wel 1 . 

As  a  minimum,  six  components  are  necessary  for  a  com- 
plete LDV  system:  laser  source,  transmitting  optics,  receiv- 
ing optics,  photodetec tor ,  signal  processor,  and  data 
processor.  The  function  of  the  transmitting  and  receiving 
optics  are  focusing  and  light  collecting,  respectively.  The 
photodetec tor  serves  to  convert  light  signal  to  an  electrical 
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signal  which  is  then  used  by  the  signal  processor  to  provide 
a  measurement  of  the  frequency.  The  capability  of  detecting 
reversed  flows  requires  an  additional  component  for  frequency 
shifting. 

Due  to  the  wide  range  of  disciplines  involved,  op- 
tics, electronics,  light  scattering,  and  signal  processing, 
the  implementation  of  this  method  can  be  difficult  CRef,  31. 
Consequently,  the  expertise  of  the  operator  can  have  a  pro- 
found impact  on  the  results  of  LDV  measurements.  It  is 
therefore  imperative,  to  an  even  greater  degree  than  many 
other  methods,  that  considerable  attention  be  given  to  ap- 
plying correct  measurement  technique. 
2.   Schlieren  FXgw  Visualization 

Schlieren  flow  visualization  systems  provide  a  quali- 
tative analysis  of  the  density  gradient  in  compressible  flow 
fields.  The  fundamental  concept  of  this  method  is  the  vari- 
ation of  the  refractive  index  gradient  of  a  fluid  with 
density . 

A  basic  schlieren  system  includes  a  light  source, 
transmitting  and  receiving  mirrors,  a  knife  edge,  and  a 
focusing  screen.  The  transmitting  mirror  serves  to  direct  a 
parallel  light  beam  through  the  test  section  and  onto  the 
receiving  mirror.  An  image  of  the  light  source  is  then 
formed  at  the  focal  point  of  this  mirror.  The  knife  edge, 
placed  in  the  plane  of  the  light  source  image,  cuts  off  part 
of  the  transmitted  light  and  controls  the  contrast  of  the 
image  focused  on  the  screen.  Changes  in  test  section  density, 
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therefore,   produce  variations  of  light  intensity  and  provide 
a  visual  representation  of  the  flow  field.  CRef.  41 

B.   INSTRUMENTATION  AND  EQUIPMENT 
1  .   Equipment  Setup 

The  experimental  procedures  for  this  research  were 
conducted  at  the  NASA  Axnes  Fluid  Mechanics  Laboratory.  A 
25-cm  by  13-cm  indraft  wind  tunnel  with  a  slotted  wall  test 
section  was  used.  Figure  2.1  shows  the  test  seption  with 
plenum  chambers  above  and  below  the  slotted  walls.  The 
adaptive  wall  design  of  this  test  section  was  developed  for 
use  in  previous  experiments.    Although  the  technique  was  not 


Figure  2.1   Wind  tunnel  test  section 
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used,  the  slotted  walls  were  helpful  in  reducing  wall 
i  nter  f erence . 

The  model  for  the  experiment  was  a  NACA  0012  airfoil 
with  3"  chord.  The  airfoil,  attached  to  the  test  section 
walls  at  quarter  chord,  was  connected  with  mechanical  linkage 
to  an  electric  motor,  as  shown  in  Figure  2.1 .  This  variable 
speed  motor  provided  a  means  for  oscillating  the  model  about 
a  pitch  axis  and  obtaining  unsteady  data  at  various 
frequencies. 

An  optical  encoder,  also  driven  by  the  motor,  pro- 
vided 3,600  counts  per  revolution  of  14-bit  data  on  the  phase 
angle  of  the  oscillating  airfoil.  The  characteristics  of  the 
encoder  were  absolute  in  that  a  whole  word  output,  with  a 
unique  code  pattern,  represented  each  of  the  3,600  discrete 
positions.  This  code  was  derived  from  independent  tracks  on 
the  encoder  disk  corresponding  to  individual  photodetec tors . 

Test  section  liach  number  and  velocity  were  obtained 
by  using  a  24-port  Scanivalve  with  a  differential  pressure 
transducer.  Calculation  of  these  values  required  meas- 
urements from  four  of  the  ports,  two  for  Mach  number  computa- 
tion and  two  for  calibration  of  the  transducer.  The  linear 
relationship  between  the  measured  pressure  and  the  corre- 
sponding voltage  produced  by  the  transducer  was  used  for 
determining  the  pressures.  In  accomplishing  this,  one  port 
of  the  Scanivalve,  vented  to  the  atmosphere,  measured  ambient 
pressure  and  a  second  port,  tapped  to  a  mechanical  vacuum 
pump,   measured  a  calibration  pressure.    A  conversion  factor 
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was  then  established  by  dividing  the  difference  of  the  known 
pressures,  atmospheric  and  calibration,  by  the  difference  in 
voltages  produced  by  the  transducer. 

Location  of  the  LDV  system  probe  volume,  or  beaim 
crossing  point,  was  controlled  by  a  three  axis  traversing 
mechanism.  A  Summit  Taskmaster  numerical  control  computer, 
rriginally  developed  for  controlling  movements  of  industrial 
machine  tools,  was  used  for  this  task.  The  numerical  control 
computer  ^ NCC  )  was  capable  of  manual  operation,  programmed 
operation,  and,  with  software  adaptation,  operation  from  the 
computer  workstation.  The  NCC  provided  control  inputs  to  two 
sets  of  electric  motors,  one  for  the  transmitting  optics  and 
the  other  for  the  receiving  optics.  The  necessity  o£  using 
different  traversing  motors  for  the  transmitting  and  receiv- 
ing optics  placed  a  great  demand  on  the  synchronization  of 
this  system.  Due  to  the  high  sensitivity  of  the  LDV  optics 
train,  the  synchronized  movement  was  often  inadequate  for 
maintaining  signal  strength. 

Phase  encoder  data  and  digitized  velocity  data  from 
the  LDV  signal  processors  were  channeled  to  a  3D-LDV  Computer 
Interface  (CI).  The  CI,  multiplexing  the  inputs  to  a  single 
digital  data  channel  output,  provided  16-bit  parallel  output 
data  to  the  computer  data  bus.  The  CI  was  capable  of  accept- 
ing up  to  six  channels  of  digital  data  and  from  one  to  three 
event  pulses  from  the  signal  processors.  The  multiplexing 
cycle  could  be  initiated  by  random  event  pulse  inputs  or 
synchronized  inputs,   the  latter  being  necessary  for  unsteady 
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data  acquisition.  When  operating  in  the  synchronized  mode, 
two  or  three  pulses,  as  selected,  occuring  within  a  set  time 
(coincidence  time)  would  be  required  to  initiate  the  multi- 
plex cycle  and  feed  data  to  the  computer. 

2.   Laser  Doppler  Velocimeter  Subsystems 

A  two  component  LDV  system  was  used  to  obtain  the 
flow  field  velocity  data.  Figures  2.2  and  2.3  show  the  gen- 
eral arrangement  of  the  system  components  for  the  transmit- 
ting  and  receiving  optics  sides,  respectively. 

A  5  Watt  Argon  ion  laser,  Spectra  Physics  Model  164, 
was   used  for   the  system.    The   unit   was   operated   in  the 


Figure  2.2   LDV  transmitting  optics 
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Figure  2.3   LDV  receiving  optics. 

multiline  mode  and  produced  a  vertically  polarized  beam  com- 
posed of  all  wavelengths  inherent  to  this  type  of  laser. 

Figure  2.4  illustrates  the  sequence  of  steps  required 
of  the  transmitting  optics.  A  set  of  mirrors  was  first  used 
to  lower  the  multiline  beam  to  the  level  of  the  optical 
components.  After  passage  through  a  collimator,  used  to 
ensure  collimation  of  the  beam  and  subsequent  crossing  of  the 
beams  at  their  waists,  the  beams  pass  through  an  attenuator. 
The  function  of  this  component  was  to  reduce  the  bean  intens- 
ity without  inhibiting  transmission  of  any  wavelengths,  a 
possible  consequence  of  simply  reducing  the  power  output  of 
the  laser.  Prior  to  separation  of  the  wavelengths  by  the 
dispersion   prism,   it  was  necessary  to  horizontally  polarize 
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Figure  2.4   Two-component  transmitting  optics 
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the  beam  by  passing  it  through  a  half-wave  plate,  or  polari- 
zation rotator,  in  order  to  obtain  maximum  light  transmis- 
sion. The  two  strons^r  w  a  veleng*:hs  ,  514.5  nm  C  green  )  and  488 
nm  (blue),  were  then  guided  by  turning  mirrors  to  the  remain- 
ing optical  c  om  p  o  n  e  n  t  s . 

Before  reaching  their  respective  Bragg  Cells,  the 
green  beam  was  split  vertically  and  the  blue  wa?'  displaced 
to  the  center  of  the  component  assembly,  vertically  polar- 
ized, and  split  horizontally.  The  relationship  between  the 
splitting  and  polarizing  planes  was  driven  by  signal  effi- 
ciency. High  signal-to-noise  ratio  requires  split  beams  to 
have  equal  intensities.  The  beams  will  be  of  approximately 
equal  intensity  when  the  polarity  of  the  incident  beam  is 
perpendicular  to  the  plane  of  the  exiting  beams  CRef.  31. 
Therefore,  the  polarity  of  the  vertically  split  beam  was 
maintained  as  horizontal  and  that  of  the  horizontally  split 
beam  was  rotated  to  vertical.  As  later  discussed,  this  dif- 
ference in  polarity  was  also  necessary  for  separation  of  the 
two  wavelengths  by  the  polarization  beamsplitter  in  the  re- 
ceiving optics. 

Bragg  Cells  are  acousto-opt ic  devices  used  to  gen- 
erate a  frequency  shift  in  the  beam  and  allow  for  resolution 
of  the  180  degree  directional  ambiguity.  One  beam  of  each  of 
the  two  wavelengths  was  directed  through  a  Bragg  Cell  pro- 
ducing a  40  MHz  shift.  The  unshifted  beams  were  passed 
through  glass  rods  to  compensate  for  the  path  length  differ- 
ence  experienced  by  the  shifted  beam.    Due  to  the  nature  of 
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the  Bragg  Cell,  the  frequency-shifted  bean  was  deflected  from 
its  initial  alignment.  Therefore,  beam  alignment  prisms,  or 
steering  wedges,  within  the  Bragg  Cell  modules  served  to 
bring  the  beams  back  to  parallel. 

A  microscope  objective,  placed  in  the  test  section  at 
the  probe  volume,  was  used  to  observe  the  beams  as  adjust- 
ments were  made  to  ensure  the  beams  focused  to  a  common 
point.  If  the  two  probe  volumes  are  not  at  the  same  location 
within  the  test  section,  then  the  measurements  will  be  made 
at  different  points.  The  necessity  for  all  four  beams  to 
overlap  at  the  focal  point  was  especially  prevalent  for  the 
acquisition  of  unsteady  data  due  to  the  need  for  a  high 
coincidence  data  rate. 

Receiving  optics  are  designed  to  image,  to  the  great- 
est extent  possible,  the  light  scattered  at  the  probe  volume 
onto  the  photodetec tors .  Since  the  scattered  light  intensity 
is  the  greatest  in  the  forward  direction,  a  forward  scatter 
system  is  most  often  preferred  for  obtaining  the  greatest 
signal-to-noise  ratio.  The  system  used  for  this  analysis  was 
an  off-axis,  forward  scatter  configuration.  Off-axis  col- 
lection angles  usually  allow  good  noise  isolation  by  admit- 
ting only  light  scattered  at  the  sample  volume  and  improve 
measurement  resolution  by  reducing  the  probe  volume  length. 
CRef.  31 

The  first  three  lenses  in  the  receiving  optics  were 
used  to  collect  scattered  light  from  the  probe  volume  and 
collimate  the  beam  (Figure  2.5).    The  beaun  was  then  directed 
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to  a  polarization  beamsplitter  for  separation  of  the  two 
colors.  The  two  components  were  then  passed  through  line 
filters,  for  blocking  ambient  light  and  preventing  cross 
talk,  and  to  lenses  used  for  focusing  the  beams  onto  the  pin- 
hole aperature  of  their  respective  photomul t i pi i er  tubes. 
The  photomul tipl ier  tubes  converted  the  light  energy  into 
electronic  signals  which  were  then  passed  to  signal  proces- 
sing equipment.  Preamplifiers  were  located  in  close  proxim- 
ity to  the  photomul t i pi ier  tubes  and  used  to  increase  the 
signal  level  before  the  introduction  of  additional  noise  into 
the  cables. 

Signal  processing  was  accomplished  with  Macrodyne 
counter  processors.  Counter  processors  operate  by  employing 
a  crystal  controlled  oscillator  to  time  over  a  selected  num- 
ber of  cycles  within  the  burst  signal.  The  clock  count  is 
initiated  when  the  filtered  signal  exceeds  a  set  threshold. 
The  Doppler  frequency  is  then  determined  by  the  time  required 
to  count  a  specified  number  of  cycles.  The  Macrodyne  proces- 
sors use  a  125  MHz  clock  and  count  over  8  or  16  cycles,  as 
se lee  ted . 

The  Macrodyne  processors  use  a  periodicity  check 
coupled  with  a  three-level  validation  scheme  to  reject  noise. 
The  periodicity  check  uses  two  counters  to  make  period  meas- 
urements over  a  different  number  of  cycles,  5  and  8  in  this 
procedure.  If  the  ratio  of  the  two  measurements  is  not  with- 
in a  specified  percentage  of  5/5,  the  signal  is  rejected. 
Three-level  validation  requires  a  signal  to  exceed  a  positive 
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Figure  2.5   Data  collection  and  processing. 
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threshold  level,  pass  through  zero,  and  then  exceed  a  nega- 
tive threshold  level  before  the  next  zero  crossing.  If  these 
conditions  are  not  met  on  any  cycle  the  counter  will  reset 
and  immediately  process  a  new  Doppler  burst.  This  method 
preven+"3  rejection  of  potentially  good  signals  having  noisy 
cycles  at  the  beginning  of  the  burst. 
3.   Schlieren  System 

The  Schlieren  system  included  two  concave  mirrors, 
each  with  a  120  inch  focal  length,  a  stroboscopic  light 
source,  a  vertical  knife  edge,  a  focusing  lens,  and  a  polar- 
oid film  support  for  obtaining  photographs  of  the  flow  field. 
The  stroboscope,  a  type  3015  Strobrite,  has  a  maximum  output 
intensity  of  300  million  beam  candles  and,  using  an  internal 
oscillator,  can  provide  flashing  rates  in  excess  of  1 ,000  per 
second.  By  adapting  a  lens  cover  with  a  vertical  slit  to  the 
stroboscope  a  slit  source  was  obtained  and  the  resolution  of 
the  image  significantly  improved. 

To  obtain  the  unsteady  photographs  a  digital  compar- 
ator was  used  to  trigger  the  stroboscope  and  provide  the 
light  for  the  photographic  imaging  on  the  polaroid  film.  The 
comparator  received  14-bit  phase  input  from  the  optical  en- 
coder and  provided  the  triggering  pulse  when  the 
instantaneous  phase  angle  of  the  airfoil  coincided  with  the 
angle  specified  to  the  comparator. 
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C.   DATA  ACQUISITION  AUTOMATION 

A  Digital  VAXstation  II  minicomputer  with  graphics  work- 
station was  used  to  automate  the  data  acquisition  process. 
The  signal  path  is  illustrated  in  Figure  2.5  starting  from 
the  probe  volume,  through  the  signal  processing  phases,  and 
terminating  within  computer  memory.  The  main  program  used  in 
accomplishing  this  task,  'LDV_ACQ',  was  originally  developed 
to  read  data  from  the  analog  to  digital  converter.  This  was 
later  modified  to  obtain  pressure  data  for  calculation  of  the 
test  section  velocity  and  call  supporting  subroutines  for  the 
acquisition  and  analysis  of  the  velocity  data. 

The   subroutines   assisting  with  various  aspects   of   the 

acquisition   process   are  listed  below  and   included   in   the 

appendix  of  this  report. 

(1  )  HWTRVT  -  allows  control  of  the  traversing  mechanism 
from  the  computer  terminal  by  sending  character  data 
to  the  interface. 

(2)  HWWA  -  combined  with  an  internally  called  subrou- 
tine, DAMUX ,  retrieves  data  from  multiplexer,  inter- 
prets phase  and  velocity  information,  creates  data 
files,  and  generates  histogram  arrays  and  scaling 
parameters . 

(  3  ■)  PT_PLOT  -  provides  graphical  display  of  computed 
velocity  values  corresponding  to  airfoil  phase 
angles  . 

(4)  PPLGT  -  similar  to  PT_PLOT  but  provides  the  addit- 
ional capability  of  allowing  the  user  to  vary  the 
plotted  velocity  range. 

(5)  HIST  -  plots  a  histogram  of  the  percentage  of  data 
points  which  fall  within  specified  velocity  limits. 
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III.   PATA  ACQUISITION 

A.   TEST  SECTION  VELOCITY 

The  main  program  for  the  acquisition  processi  •LDV_ACQ', 
was  responsible  for  the  calculation  of  test  section  velocity. 
In  combination  with  scanivalve  positioning  subroutines,  the 
program  acquired  pressure  transducer  data  from  7  ports.  As 
previously  discussed,  measurements  from  2  of  the  ports  were 
used  for  calibration.  Pressure  data  from  3  ports  was  dis- 
carded, and  the  other  two  provided  the  static  pressure  at  the 
entrance  to  the  test  section  and  the  total  pressure  measured 
at  the  nozzle  inlet.  Assuming  isentropic,  perfect  gas  condi- 
tions, this  supplied  the  necessary  parameters  for  the  Mach 
number  computation  by  employing  the  following  relationship: 

M=(((P,/P)^'-'''-l)*(2/(7-l)))'^' 
with  a  value  of  1 .4  used  for  7  •    The  test  section   velocity 

was  then  calculated  by: 

V  =  H  X  a 

where  a,  the  sonic  velocity,  was  obtained  as  follows: 


a  =  VyRT 
The  temperature  required  for  this  calculation  was  read  from  a 
wall   thermometer  in  the  test  cell  and  entered  when   prompted 
by  the  program. 

The   velocity   values  were  transferred  to  data  files   and 
later  used  to  normalize  the  LDV  velocity  data. 
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B.  SCHLIEREN  PHOTOGRAPHY 

A  qualitative  analysis  of  the  steady  and  unsteady  flow 
fields  was  conducted  at  two  Mach  numbers,  0.3  and  0.5,  and 
angles  of  incidence  varying  from  0  to  20  degrees.  In  ac- 
complishing this,  test  section  Mach  number  was  set  by  adjust- 
ing the  wind  tunnel  throat  and  verified  with  Scanivalve 
pressures  and  calculations  from  'LDV_ACQ'.  Final  adjustments 
to  the  Stoboscope  intensity  and  knife  edge  were  then  made  to 
ensure  optimum  contrast. 

C.  WAKE  VELOCITIES 

With  the  desired  wind  tunnel  conditions  established  and 
the  probe  volume  in  position  for  the  data  point,  'HWWA',  a 
subroutine  called  by  the  main  program,  was  employed  for  raw 
data  acquisition  and  conversion  to  velocity  and  phase  infor- 
mation. An  internally  called  subroutine,  'DAMUX',  acquired  4 
or  6  channels  of  data,  as  selected,  from  the  LDV  Computer 
Interface  (CI).  The  CI  was  capable  of  transferring  either  2, 
4,  6,  or  8  channels  of  data,  with  the  first  two  of  the 
channels  reserved  for  time  and  status  information  generated 
by  the  CI.  During  the  acquisition  of  steady  state  velocities, 
4  channels  were  selected  to  provide  the  2  velocity  components 
in  addition  to  time  and  status  data,  which  were  discarded 
since  they  were  not  relevant  for  the  experiment.  For  unsteady 
acquisition,  6  channels  of  data  were  transferred  to  accom- 
modate the  additional  data  channel  required  for  phase   input. 
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As  previously  mentioned,  the  signal  processors  determined 
the  Doppler  frequency  by  using  the  time  required  to  count  a 
specified  number  of  cyles,  8  for  this  analysis.  This  raw 
data  was  transferred  to  the  computer  in  16  bits,  with  a  10 
bit  mantissa,  a  4  bit  exponent,  and  2  nonessential  bits  used 
for  controlling  data  flow.  Due  to  the  format  of  the  time 
words,  it  was  necessary  to  invert  and  reverse  the  bits  before 
extracting  the  exponent  and  mantissa  values.  The  signal 
period  was  then  calculated  as  follows: 

ts  =  < 1/32)< D9D8.  .  .DO  )  x  2 
where   D9D8...D0  represented  the  binary  mantissa  and  EXP   was 
the   range   selected   with  the  processor,  used  to  divide   the 
number  of  clock  counts  in  the  Doppler  burst. 

The  signal  frequency,  the  inverse  of  the  period,  was  ob- 
tained from: 

f,  =  (32  X  1000)/<2^^^  X  mantissa) 

The  signal  frequency  was  composed  of  the  Doppler  frequen- 
cy, the  40  MHz  shift  frequency,  and  the  downmixing  frequency. 
This  allowed  the  Doppler  frequency  to  be  deduced  as  follows: 

4  t     thifl    ml< 

The  8  cycles  counted  within  the  Doppler  signal  physically 
represent  8  fringes  crossed  by  a  particle  travelling  with  the 
flow.  Consequently,  the  velocity  component  perpendicular  to 
the  fringe  pattern  may  be  calculated  by  the  product  of  the 
Doppler  frequency  and  the  fringe  spacing.  The  equation  for 
the  latter  is: 
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<5  =  X  /  [2  *  sin(d/2)] 
where  A   is  the  wavelength  of  the  beam  and  9    is  the  angle  be- 
tween the  intersecting  beouns . 

After  obtaining  a  sufficient  number  of  data  points  at  a 
probe  position  (2,000  for  unsteady  flows),  the  velocity  and 
phase  data,  as  well  as  the  probe  position,  test  section 
velocity,  and  mixing  frequencies,  were  stored  in  files  for 
later  analysis.    The  Taskmaster  traverse  controller  was  then 

■ 

used  to  move  the  probe  volume  and  data  collection  continued. 

The  data  acquisition  software  included  two  routines  for 
on-line  performance  checks  of  the  data.  'HIST*  was  used  to 
provide  a  histogram  of  the  percentage  of  data  points  which 
occurred  within  specified  velocity  limits.  Figure  3.1  shows 
a  sample  plot  provided  by  the  subroutine.  This  data  valida- 
tion method  is  most  useful  for  steady  flow  conditions,  in 
which  case  a  high  percentage  of  data  concentrated  around  the 
estimated  velocity  would  be  a  positive  indicator.  For  un- 
steady analysis  'PT_PLOT'  provided  a  graph  of  computed  ve- 
locity values  corresponding  to  airfoil  phase  angles.  With 
this  information  available,  the  data  quality  could  be 
easily  assessed  and  a  determination  made  regarding  the  ac- 
ceptance of  these  values. 
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Figure  3.1   Sample  histogram 
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IV.   RESVLTg  Am   DISCUSSION 

A.   SCHLIEREN  PHOTOGRAPHS 

The  Schlieren  photographs  obtained  from  the  analysis  of 
the  steady  flow  field  were  intended  to  provide  information  on 
the  parameters  corresponding  to  early  flow  separation  and  a 
reference  for  comparison  with  unsteady  flows  at,  the  same 
conditions.  Steady  state  separation  was  first  observed  at  a 
Mach  number  of  0.5  and  a  10  degree  angle  of  incidence  (Fig- 
ure  4.1  ).    Based  on  this,   the  parameters  for  the   unsteady 


Figure  4.1   Steady  state  Schlieren  photograph; 
mean  a  =1 0  degrees. 
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analysis  were  selected  to  be  a  Mach  number  of  0.5  and  a  phase 
angle  corresponding  to  a  1 0  degree  angle  of  incidence.  Un- 
steady flow  photographs  were  then  obtained  at  oscillation 
frequencies  of  5  and  25  Hz  for  comparison  with  each  other  and 
with  wake  profiles  generated  from  LDV  data  (Figures  4.2  and 
4.3). 

From  a  comparison  of  the  steady  and  unsteady  photographsi 
it  would  appear  that  oscillating  the  airfoil  reduces  the  size 
of  the  separated  wake  and,  consequently,  mitigates  the  effect 
of  stalling.    This  is  more  apparent  at  the   higher  frequency 


Figure  4.2   Unsteady  flow  Schlieren  photograph; 
f=5  Hz,  a  =10  degrees. 
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Figure  4.3   Unsteady  flow  Schlieren  photograph; 
f=25  Hz,  a  =10  degrees. 


oscillation,   a   logical  result  due  to  the  greater  effect   of 
time  unsteadiness  at  higher  reduced  frequencies. 

B.   DATA  ANALYSIS  SOFTWARE 

The  main  program,  'ORCIZE*,  and  its  supporting  subrou- 
tines were  designed  to  provide  wake  profile  plots  from  the 
data  files  generated  during  acquisition.  To  accomplish  this, 
the  data  file  corresponding  to  each  point  in  the  profile  is 
input  to  the  software.  A  series  of  routines  first  reorgan- 
izes the  data  files  and  discards  statistically  bad  data 
points,   those   outside   of  two  standard  deviations  from   the 
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mean.  Due  to  the  oscillatory  nature  of  the  forcing  function 
producing  the  unsteady  data,  a  sinusoidal  curve  fit  is  then 
introduced  to  allow  evaluation  of  the  first  harmonic  of  the 
velocity  distribution  over  the  phase.  A  unique  sinusoid  is 
generated   for  each  data  point  within   the   profile. 

A  minor  degree  of  uncertainty  is  introduced  into  the  ve- 
locities by  the  data  interpolation  and  curve  fitting  rou- 
tines discussed  above.  However,  unlike  many  other  velocity 
measurement  techniques,  LDV  data  is  obtained  randomly  due  to 
the  random  arrival  of  the  seed  particles  at  the  prove  volume. 
Without  this  curve  fit,  it  would  be  difficult  to  generate 
velocity  information  at  the  same  desired  phase  angle  for  each 
data  point  in  the  profile. 

During  the  generation  of  the  wake  profile  the  user  is 
provided  a  comparison  plot  of  mean  values  and  the  sinusoid 
curve  fit.  The  mean  values  are  provided  by  'BDPTS'  and  com- 
puted within  each  20  degree  phase  window,  as  illustrated  in 
Figure  4.4.  This  was  helpful  in  ascertaining  the  quality  of 
the  curve  fit  and,  consequently,  the  quality  of  the  wake  pro- 
file results. 

The  data  analysis  routines  are  summarized  below  and  in- 
cluded in  the  appendix  of  this  report. 

(1  )  ORCIZE  -  main  progrsim ,  calls  supporting  subroutines; 
sequentially  organizes  data  to  assist  with  subse- 
quent analysis. 

(2)  PT_PLOT  -  provides  graph  of  computed  velocity  values 
corresponding  to  phase. 

(3)  BDPTS  -  computes  mean  and  standard  deviation  of  ve- 
locities  occurring  within  20  degree  phase   windows; 
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Figure  4.4   Comparison  of  mean  and  curve  fit 


discards  data  points  outside  of  two  standard  devia- 
tions of  the  mean  within  each  of  these  windows;  pro- 
vides standard  deviation  plots  of  data. 

(4)  AVG  -  averages  data  values  occurring  at  same  dis- 
crete phase. 

<5)  INTERP  -  provides  a  linear  interpolation  between 
data*  resulting  in  180  equally  spaced  points;  neces- 
sary for  input  to  FOURIER. 

(6)  FOURIER  -  computes  coefficients  for  a  sinusoidal 
curve  fit  of  the  data . ^ adapted  from  June  1977  OCS 
Airfoil  program). 

(7)  CVFIT  -  fits  sinusoidal  curve  to  velocity  data  ( u- 
tilizes  FOURIER). 


34 


(8)  PLOT   -  provides   graphical  comparison  of  curve   fit 
with  mean  values  generated  by  BDPTS . 

(9)  PROPLOT  -  provides  wake  profile  plots. 


C.   WAKE  PROFILES 

LDV  data  was  obtained  for  steady  and  unsteady  flows  at  a 
riach  number  of  0.5  and  a  10  degree  angle  of  incidence.  The 
steady  profiles  were  collected  at  distances  of  0.5,  1.0,  and 
2.0  inches  behind  the  trailing  edge  of  the  airfoil.  These 
values  were  non-dimensional ized  by  the  chord  length  of  the 
airfoil  and  are  indicated  as  such  on  the  profile  plots.  Due 
to  time  limitations,  only  one  data  profile  was  collected  for 
the  unsteady  case,  at  two  inches  behind  the  trailing  edge. 

The  wake  profile  plots  are  presented  in  the  following 
pages  of  this  section.  A  lack  of  data  points  in  the  lower 
part  of  the  flow  field  will  be  apparent  upon  inspection  of 
these  results.  This  was  due  to  a  mechanical  restriction  on 
the  traversing  mechanism. 

The  curve  fit  of  the  plotted  values  of  the  steady  data 
was  approximated  and  hand  drawn  to  aid  the  reader  in  ob- 
serving the  effects  of  the  airfoil  in  the  flow  field.  The 
software  provides  a  spline  interpolation  fit  between  the 
data,  as  used  in  the  unsteady  plots.  However,  due  to  the 
amount  of  data  scatter  involved,  this  method  was  not  used  for 
these  results. 
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HORIZONTAL  VELOCITY  PROFILE 
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Figure  4.5   Steady  wake  profile 
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VERTICAL  VELOCITY  PROFILE 


on 


CJ> 


cc: 


CO 


CD 

en 


e.  .  uu    - 

X/C  =  0.2 

i 
• 

\ 

1.50   - 

• 
1 

• 
/ 

1   00  - 

• 

\ 

• 

1 

• 

0.50   - 

• 
• 

0   00  - 

^^      • 

-0.50  - 

-1    00  - 

-1.50   - 
-2.00  - 

r               r 

1 ,■              , 

-0.70    -0.50    -0.30    -0.10     0.10     0.30     0.50     0.7i 

VELOCITY  RATIO  (V/U(INFIN)) 


Figure  4.6   Steady  wake  profile 
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Figure  4.7   Steady  wake  profile, 
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Figure  4.8   Steady  wake  profile 
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Figure  4.9   Steady  wake  profile 
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Figure  4.10   Steady  wake  profile, 
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Figure  4.11   Unsteady  wake  profile 
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Figure  4.12   Unsteady  wake  profile 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  purpose  of  the  results  obtained  during  the  experi- 
mental phase  of  this  thesis  was  to  validate  the  technique  and 
software  developed.  This  was  accomplished  when  unsteady 
data,  generated  by  an  oscillating  airfoil,  was  successfully 
acquired,  analyzed,  and  reduced  to  wake  profile  plots. 

Throughout  the  course  of  this  investigation  hardware 
problems  have  a  posed  a  signigicant  challenge  to  acquiring 
data  within  the  existing  time  constraints.  Consequently,  the 
data  obtained  was  accepted  and  included  in  this  report  with- 
out the  opportunity  to  more  carefully  investigate  regions  of 
significant  data  scatter  and  improve,  or  validate,  the  data 
plots . 

The  results  produced  by  these  procedures  confirm  a  sound 
technique,  but  also  reveal  several  areas  in  need  of  addi- 
tional emphasis.  At  present,  an  interpolation  subroutine  is 
required  to  produce  data  at  equally  spaced  intervals  before 
generation  of  the  curve  fitting  coefficients  by  'FOURIER*. 
Combined  with  the  inability  to  utilize  a  large  number  of  data 
points  (  <  300),  this  limitation  will  inherently  introduce 
additional  uncertainty  into  the  results.  A  second  software 
modification  to  be  considered  is  the  generation  of  the  veloc- 
ity data  at  various  phase  angles  simultaneously.  Currently, 
the   data   must   be  processed  for  each  separate   phase   angle 
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desired.  This  would  eliminate  the  need  for  storage  of  large 
amounts  of  raw  data  and  significantly  reduce  the  processing 
time.  By  increasing  the  range  of  vertical  movement  of  the 
traverse,  the  mechanical  restriction  that  prevented  the  ac- 
quisition of  data  in  the  lower  regions  of  the  wake  can  be 
el iminated . 

Continuation  of  this  study  will  be  accomplished  with  a 
newly  developed  test  section  currently  being  installed.  This 
test  section  will  enable  schlieren  flow  visualization  and  LDV 
data  acquisition  at  locations  over  the  airfoil  which  were 
previously  blocked  by  the  mechanical  linkage  used  for  sup- 
porting and  oscillating  the  airfoil.  This,  combined  with  the 
modifications  suggested  above,  will  enable  an  in-depth  analy- 
sis of  wake  and  boundary  layer  characteristics  associated 
with  the  dynaimic  stall  phenomenon. 
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APPENDIX 
COMPUTER    PROGRAMS 


PROGRAM  LDV  ACQ 


•  Adapted  froB  LABSTAR  Exaaple3_l  by  S.  D.  Hedrick  in  Hay,  1987  for  reading   • 

•  A-0  data.   Modified  in  August,  1987  to  be  used  as  the  central  progcaa  in    • 

•  the  acquisition  of  unsteady  velocity  coaponents  with  an  LDV.  • 

•  • 

INCLUDE  'sysSlibrary:LIOSET.rOR' 
C  declare  local  variables 

INTEGER  sys_stat         ! Status  ret  by  LIO  calls 

INTEGER  dev_id,  clk_id   ILIO  device  ID 

INTEGER  data_length      Inuaber  of  bytes  of  data  read 

INTEGER  i,  npt,  dnpt 

INTEGER  chn,  nsteps,  ier 

REAL     rate,  volt(0:23),  press(0:23),  corfac 

REAL    sua,  aean,  calp,  conv,  aab,  teap,  aach,  velo 

CHARACTER* 1  ans 

CHARACTER*4  clock 

LOGICAL  hoae 

C  declare  data  buffers 

INTEGER*2  raw_data ( 501 ) 
REAL    VOLTAGES ( 501  ) 


C 


OPEN (UNIT- 13, FILE-' scani .dat' , STATUS- 'NEW' ) 

write(*,*)  'ENTER  A-D  CHANNEL  TO  BE  USED  (0,1,8,9)' 
read( • , *  )  chn 
10     write(*,*)  'ENTER  CLOCK  RATE  (.01  -  4500  H2)--ALL0W  FOR  A', 
1  '  MAXIMUM  OF  2  SECONDS  TO  ACQUIRE  DESIRED  SAMPLES' 

read(*,*)  rate 
if  (rate  .It.  .01  .or.  rate  .gt.  4500)  then 

write(*,*)  'RATE  IS  NOT  WITHIN  SPECIFIED  LIMITS' 
go  to  10 
end  i  f 
20     write(*,*) 'ENTER  NUMBER  OF  DATA  POINTS  TO  BE  COLLECTED', 
1  '  (500  aax) ' 

read(*,*)  npt 

if  (npt  .It.  1  .or.  npt  .gt.  500)  then 

write(*,*)  'NUMBER  OR  POINTS  IS  NOT  WITHIN  LIMITS' 
go  to  20 
end  if 

npt  -  npt  ♦  1 
dnpt  -  npt  •  2 

writeC,*)  'ENTER  CLOCK  DEVICE  (KZAO.KZBO)' 
read(*, ' (A4)' )  clock 
C 
1       C  Attach  the  ADVll-D  and  set  up  for  queued  I/O 
!       C 

I  sys_Stat-LIOSATTACH(dev_id, 'AZAO' ,  LIO$K  QIO) 

I  IF(.NOT.  sys  Stat)  CALL  LIBSSIGNAL( %val ( sys_stat ) ) 

I       C 

C  Set  up  the  ADV  for  synchronous  transfer  on  selected  channel 

C        Gain  of  1 

C 

sysstat    -    LIOSSETI    (dev    id,    L10$K    SYNCH,    0) 

IF(.NOT.    sys    statT   CALL    LIB?SIGNAL( %val ( sys    stat)) 
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sys_6tat  -  LIO$SET_I  ( dev  id,  LIO$K  AD_CHAN,  1,  chn) 

IF( .NOT.  sys_statT  CALL  LIB$SIGNAL( %val ( sys_stat ) ) 

sys_stat  -  LIO$SET_I  (dev  id,  LIOSK  ADGAIN,  1,  1) 

IFI.NOT.  sys_statT  CALL  LIB$S1GNAL( *val ( Bys_6tat ) ) 
C 

C  Get  a  rawdata  buffer 

C       length  of  buffer  and  returned  data  length  is  in  bytes 
C  Attach  clock  and  set  rate 
C 

hose  -  .true. 

nsteps  -  0 

do  40  j-0,7 

call  dalljllab(hoBe , nsteps, ier ) 

sys_stat-  LIO$ATTACH(clk_ id, clock, LIO$K_QIO) 

ir( .NOT.  sys  Stat)  CALL  LIB$SIGNAL( Ival ( sys_stat ) ) 
sys_stat-LIO$SET_R(cIk_id,LIO$K_CLK  RATE,  l,rate) 

IF( .NOT.  sys  Stat)  CALL  LIB?S1GNAL( %val { sys_stat ) ) 
sys_Stat-LIO$SET_I(cIk_id,LIO$K_TRIG,  1,  LIO$K  IMMEDIATE) 

ir( .NOT.  sys  Stat)  CALL  LIB$SIGNAL( %val( sys_stat ) ) 
sys_stat  -  LI0$REA0(3ev  id,  raw_data ,dnpt ,data  length,  ) 

ir( .NOT.  sys_stit)  CALL  LIB$SIGNAL( %vaT( sys  stat)) 
sys  Stat  -  LIOSSET  I  (dev_id,  LIO$K_TRIG,  2,  LIO$K  CLK  POINT, 
1  clk_idT 

IP( .NOT.  sys_stat)  CALL  LIB$SIGNAL( %val ( sys  Stat)) 
C 

C  Convert  the  raw  data  to  voltages 
C 

call  LSP$FORMAT_TRANSLATE_ADC( raw_data, voltages, npt,  ,  ,  ) 
C 
C 

C  calculate  the  mean  value  of  the  voltages  and  detach  the  clock 
C 

sum  >  0.0 

do  30  i-2,npt 

sum  «  sum  -f  voltages(i) 
30     continue 

mean  -  sum  /  float(npt-l) 

vol t ( j )  •  mean 

sys_stat  -  LIO$DETACH(clk  id,  ) 

IF(.NOT.  sys_statT  CALL  LIB$S1GNAL( %val ( sys_stat ) ) 
40     continue 
C 

C  detach  from  the  A/D 
C 

sys_stat  -  LIO$DETACH(dev  id,  ) 

ir(.NOT.  sys  statT  CALL  LIB$SIGNAL( Ival ( sy6_stat ) ) 
C 

C  convert  scanivalve  voltages  into  units  of  pressure 
C 

wrlte(«,*)  'ENTER  CALIBRATION  PRESSURE  (in.  of  Hg)' 
read(*,*)  calp 

write(*,*)  'ENTER  REQUIRED  CORRECTION  FACTOR  FOR  CALIBRATION', 
1  '  PRESSURE  (in.  of  Hg ) ' 

r*ad(*,*)  corfac 
calp>calp-fcorfac 

wrlt*(*,*)  'ENTER  ATMOSPHERIC  PRESSURE  (in.  of  Hg ) ' 
read(*,*)  aab 
prtst(O)  -  calp 
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press( 1 )  ~  amb 

conv  -  (pres6(l)-pres6( 0)  )  /  ( volt ( 1 ) -volt ( 0 ) ) 

do  50  j-2,7 

pres6(j)  -  conv  •  volt(j)  ••■  calp  -  volt(O)  •  conv 
50     continue 

WRITE!  •, * )  VOLT( 0) ,V0LT(1 ) ,V0LT(2) ,V0LT(7) 

write( 13,55) 
55     for«at( 2X, 'PORT' ,4X, 'PRESSURE' ,//) 

do  60  i-0,7 

write(13, '(3X,l2,5X,r9.4)' )  i,  press(i) 
60     continue 
C 

C  coapute  test  section  velocity 
C 

write! *,*) 'ENTER  AMBIENT  TEMPERATURE  (deg.  F)' 

read! • , * )  teap 

Bach-sqrt(abs(5.*( ( press( 2 )/press( 7 ) )••( .4/1.4)-!. ) ) ) 

velo-Bach*5qrt(1.4»53.3*32.174*(teBp*4  59.67) ) 

write! *,*) 'TEST  SECTION  HACH  NUMBER  AND  VELOCITY — ',aach,velo 
C 

write! *,*) 'CONTINUE  WITH  LDV  DATA  ACQUISITION?  (Y/N,  DEF-Y)' 

read! •, ' !A1) ' )  ans 

if  (ans  .eq.  'n'  .or.  ans  .eq.  'N' )  go  to  70 
C 

C  call  subroutine  to  aove  LDV  traversing  aechanisa  and  acquire  data 
C 

call  HWTRVT(velo) 
C 

70     stop 

end 
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SUBROUTINE  HWWA( XPOS , ZPOS , VELO ) 
C 
C  ACRONYM:   Hardware  checkout  -  test  of  DRVll-WA 

C  PURPOSE: 

C        For  LVIS,  HWDVYB  drives  the  acquisition  of  the  LDV  data. 

C        It  is  a  test  program  for  the  YB  driver. 

C 

C  HETHOD: 

C 

C   BEGIN 

C   Install  loadable  drivers  if  necessary. 

C   Read  the  various  instrument  types  in  succession,  skipping  those  which 

C      are  not  to  be  read. 

C   Reaove  loadable  drivers  if  necessary. 

C   END 

C 

C  PARAMETERS: 

C     ARC        DIM       TYPE  I/O/S  DESCRIPTION 

C 

C  COMMONS  USED: 

C 

C  FILES  USED:   None 

C 

C  ERROR  HANDLING: 

C     lER  -  1  means  no  "fatal"  error. 

C 

C  NOTES: 

C 

C  LOCAL  VARIABLES: 

C     VAR        DIM       TYPE     DESCRIPTION 

C 

C  EXTERNAL  REFERENCES: 

C     NAME       DESCRIPTION  AND  SOURCE 

C     DVWA       Calls  DRVll-WA  device  driver  to  acquire  A/D  data. 

C 

C  STANDARDS  VIOLATIONS:  None. 

C 

C  ENVIRONMENT:   DEC  VAX/VMS  AND  FORTRAN  77 

C 

C  DEVELOPMENT  HISTORY: 

C        DATE    INITIALS    DESCRIPTION 

C      07/26/83    TML      Adapted  from  DAPNT 

C      12/12/85    CLH     Modified  to  print  results  in  either 

C  octal  or  integer. 

C     02/25/87    GBG     Modified  for  use  with  DRVll-WA  and  VMS 

C      05/06/87    CLH      Modified  to  accept  "new"  DAMUXl . 

C     06/08/87    RRR     Modified  to  subroutine  ♦  deleted  CLH  5/87  Mod 

C   JUL-AUG,1987   SDH      Modified  to  calculate  velocities  using  variable 

C  mixing  frequency  as  user  input,  interpret  phase 

C  angle  information  generated  by  digital  encoder, 

C  transfer  test  section  velocity  and  probe  posit- 

C  ion  to  data  file,  and  generate  histogram  arrays 

C  and  scaling  parameters 

C 

C  AUTHOR(S):  Ted  Lichtenstein,  Informatics  General  Corp. 

C 

C- •**•*•••••••••••*•••••••••••••••••••••••**•****•••*••••••*•••••*•**•*• 

C 

IMPLICIT  NONE 
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100 


200 


PARA/1ETER  NRAWHX-80000 

PARAMETER  LUNTI-5 

PARAMETER  LUNOUT-6 

INTEGER*2  I RAWDT ( NRAWMX ) ,  I fORM ,  MUXTIME  (16384) 

INTEGER*2  MUXSTS  ( 16 384 ), RAWDA( 20000 ), RAWDB( 20000  ) 

INTEGER«2  RAWDC( 20000) 

INTEGER  STATUS, I ER.PHAN 

INTEGER  NCOLUn,  lYESNO,  NRAWIN.I.J 

INTEGER  LIBSSUBX,  NCHANS ,  NSAMPLES,  INDEX,  JNDEX 

INTEGER  nRAWA,HRAWB,NRAWC,nRAWD 

I NTEGER  I RAWDTA , I RAWDTB , I RAWDTC , I RAWDTD 

INTEGER  IRAWDTE.IRAWDTF 

REAL  FREOC  2) ,VEL( 2) , HIXA, MIXB , XPOS , ZPOS , VELO 

LOGICAL  DATSELCT  (3) 

CHARACTER*10  INFILE 

CHARACTER*1  ANS 

REAL* 4  YU(160) ,yV(120) , YUK 160 ) , YVl ( 120 ) ,YMAXU,YMAXV 

REAL* 4  SAVE 

INTEGER  COUNT,XLOW,XHI 

REAL  X(4000),Y(4000) 

INCLUDE  ' (SSYSSRVNAH) ' 

WRITE( *,*) 'ENTER  DESIRED  NAME  Of  DATA  TILE — E.G.,  LDVXXX . DAT , ' , 
1  '  WHERE  THE  'XXX*  HAY  BE  USED  AS  A  SERIAL  ID  FOR  THE', 

1  '  PROBE  POSITION' 

READ( *, ' (AlO) ' )  INFILE 

open  ( uni t-11 , file-INFILE, status"' new' ) 

CONTINUE 

Fill  the  raw  data  buffer  with  this  bit  pattern: 

"1010101010101010" 

This  equals  125252(8),  -21846(10). 

DO  2  00  I  -  1,  NRAWMX 

IRAWDT(I)  -  '125252'0 
CONTINUE 


DO  I  -  1 ,  16384 


END  DO 


MUXTIME  (I)  -  •125252'0 
MUXSTS  (I)  -  '125252'0 


220  CONTINUE 

*   Prompt  for  no.  of  words  to  acquire: 

WRITE(*,*)'  ' 

WRITEC,  •) 'Starting  LDV  acquisition' 

WR1TE(*,*)'  ' 

WRITE(*,* ) 'ENTER  NUMBER  OF  CHANNELS  TO  BE  USED', 
1  '  (assuaes  non-oscillating  airfoil  if, 

1  '  lets  than  six)' 

REAO(*,*)  NCHANS 
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224  WRITE  (♦,  225) 

225  rORNAT  ('  Enter  number  of  samples:') 
READ  (*,  »)  NSAWPLES 

IF  (NSAHPLES  .GT.  20000)  THEN 
WRITE  (*,  226) 

226  rORHAT  ('  The  number  of  samples  must  be  20,000  or  less.') 

GO  TO  22  4 
END  ir 

C 

2  30   CONTINUE 
C   *   Prompt  for  output  format  of  data; 
C 

IFORM-2 

IF  (  IFORM.LT.l  .OR.  I  FORM . GT . 2  )  GO  TO  2  30 
C 
C 
C 

C       Select  the  data  to  be  returned. 
C 

DATSELCT  (1)  -  .false. 
DATSELCT  (2)  -  .false. 
DATSELCT  (3)  -  .TRUE. 
C  - 

C       Prompt  for  the  mixing  frequencies  selected  by  user 
C 

WRITE( *,*)  'ENTER  MIXING  FREQUENCIES  FOR  EACH  COMPONENT  -  (blue, 
Igreen)  -  in  MHz' 
READ(«,»)  MIXA.MIXB 

C 

C       Read  and  display  the  data  one  sample  at  a  time. 

C 


C 


IF  (NCHANS  .LT.  6)  THEN 

WRITE( *,») 'ENTER  AIRFOIL  PHASE  ANGLE  (tenths  of  degrees,  14)' 
READ(»,*)  PHAN 
DO  9  I-1,NSAMPLES 

CALL  DAMUX ( NCHANS , 1 , DATSELCT , MUXTIME , MUXSTS , IRAWDT , I ER ) 
RAWDA( I  )-IRAWDT(  1) 
RAWDB( I )-IRAWDT(2) 
RAWDC( I )-PHAN 
9     CONTINUE 
GO  TO  15 
END  IF 

DO  10  I-1,NSAMPLES 

CALL  DAMUX  (NCHANS,  1,  DATSELCT,  MUXTIME,  MUXSTS, 
•  IRAWDT,  lER) 

RAWDA( I )-IRAWDT( 1 ) 
RAWDB( I )-IRAWDT( 2) 
RAWDC( I )-lRAWDT( 3) 
10  CONTINUE 

DO  19  I-1,NSAMPLES 
HRAMA-O 
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HRAWB-0 
MRAWC-0 
MRAWD-0 

MRAWA-IIBITS(RAWDC( I ) ,12,2) 
MRAWB- 1 1 B I TS ( RAWDC ( I ) , 8 , 4 ) 
HRAWC-IIBITS{RAWDC(I ) ,4,4) 
MRAWD-IIBITS(RAWDC(I ) ,0,4) 

RAWDC( 1 )-HRAWA*1000+nRAWB*100+HRAWC*10*MRAWD 
19  CONTINUE 

15  DO  20  I-1,NSA«PLES 

RAWDAI I)-INOT(RAWDA( I) ) 
RAWDB(I)-INOT(RAWDB( I) ) 
IRAWDTA-0 
IRAWDTB-0 
IRAWDTC-0 
IRAWDTD-0 
IRAWDTE-0 
IRAWDTF-O 

IRAWDTA-IIBITS(RAWDA( I ) ,10,4) 
IRAWDTB-IIBITS(RAWDA( I ) ,0,10) 
I RAWDTC- 2  * • 1 RAWDTA • I RAWDTB 
ir    (IRAWDTC  .EQ.  0)  THEN 
WRITE( *,* ) ' IRAWDTC-0' 
GO  TO  20 
END  ir 

FREOd  )-((32.*10.**3.  )/(rLOAT(  IRAWDTC)  )  )-40.-»-MIXA 
IRAWDTD-IIBITS(RAWDB( I ) ,10,4 ) 
IRAWDTE-IIBITS(RAWDB( I ) ,0,10) 
'  IRAWDTF-2**IRAWDTD*IRAWDTE 
IF  (IRAWDTF  .EQ.  0)  THEN 
WRITE! *, •) ' IRAWDTF-0' 
GO  TO  20 
END  IF 

FREQ(2)-( (32.*10.**3.)/( FLOAT ( IRAWDTF) ) )-4  0.+MIXB 
vel(l)-freq(2)*0.514  5/(2.*sin(atan(13./(2.*482.6)))) 
vel(2)-f req( 1 ) 'O . 4880/( 2 . 'sirK  atan( 13./(2.*482.6)))) 
veil l)-vel(l)/.3048 
vel(2)-vel(2)/.3048 

WRITE (11, 7 00)  vel(2) ,vel( 1 ) , RAWDC ( I ) 
20  CONTINUE 

WRITE( 11 ,900)  XPOS,ZPOS 
WRITE( 11 ,950)  VELO,NSAMPLES 
WRITE( 11,975)  MIXA,niXB 

WRITE(*,* ) 'DO  YOU  DESIRE  A  HISTOGRAM  OF  THE  VELOCITIES  OBTAINED?' 
1  '  ( Y/N  -  DEF-N) ' 

READ( *, ' (Al ) ' )  ANS 
IF  (ANS  .NE.  'y'  .AND.  ANS  .NE.  'Y')  THEN 

GO  TO  3  5 
END  IF 

CLOSE(UNIT-ll) 

0PEN(UNIT-11,FILE-INFILE, STATUS- 'OLD' ) 
DO    40    I-1,NSAHPLES 

READ(11,800)    X(I),Y(I) 
40      CONTINUE 
XLOH--200 
XHI  — 195 
DO   60    J-1,160 
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COUNT- 0 

DO  50  I-l,NSArtPLES 

IF  (X(I)  .GE.  XLOW  .AND.  X(I)  . LT .  XHI )  THEN 

COUNT- COUNT* 1 
END  IF 
50     CONTINUE 

YU( J )-( rLOAT( COUNT )/NSA«PLES) •100. 
XLOW-XLOW-fS 
XHI-XHI-fS 
60   CONTINUE 
C 

C  Calculating  a  scaling  parameter  for  the  horizontal  velocity  histograi 
C 

DO  62  1-1,160 
YUK  I  )-yU(I) 
6  2   CONTINUE 

DO  63  1-1,160 
DO  65  J-1,159 

ir    (YU1(J*1)     .LT.     YUK  J))    THEN 
SAVE-YUK  J) 
YUK  J)-YU1(  J*l) 
YUK  J  +  1  )-SAVE 
END    ir 
65  CONTINUE 

63       CONTINUE 

YMAXU-YUK  160) +  10.0 
C 

XLOW--300 
XHI--295 
DO  80  J-1,120 
COUNT- 0 
DO  70  I-1,NSAMPLES 

IF  (Yd)  .GE.  XLOW  .AND.  Yd)  .LT.  XHI)  THEN 

COUNT-COUNT+1 
END  IF 
70     CONTINUE 

YV{ J)-( FLOAT(COUNT)/NSAMPLES)*100. 
XLOW-XLOW+5 
XHI-XHlf5 
80   CONTINUE 
C 

C  Calculating  a  scaling  parameter  for  the  vertical  velocity  histogram 
C 

DO  82  1-1,120 
YVK  I)-YV(  I  ) 

82  CONTINUE 

DO  83  1-1,120 
DO  85  J-1,119 

IF  (YVKJ^l)  .LT.  YVKJ))  THEN 
SAVE-YVK  J) 
YVK  J)-yVl(  J  +  1) 
YVK  Jt-1) -SAVE 
END  IP 
8  5     CONTINUE 

83  CONTINUE 

YHAXV- YVK  1 20  )  •4-1 0  .  0 
C 

CALL  HIST(YU,YV,yMAXU,YMAXV) 
C 

35   CLOSE(UNIT-ll) 
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WRITE(*,* ) 'DO  YOU  DESIRE  A  PLOT  Of  THE  COMPUTED  VELOCITY', 
1  •  VALUES?  (Y/N  -  DEF-Y)' 

READ( •, • (Al ) ' )  ANS 
IF  (ANS  .EQ,  'N'  .OR.  ANS  . EQ .  'n')  THEN 

GO  TO  30 
END  If 

CALL  PT_PL0T(NSAJ1PLES,INfILE) 

WRITE( *, • ) 'DO  YOU  DESIRE  TO  VARY  THE  PLOTTED  VELOCITY  RANGE?', 
1  '   (Y/N — DEF-Y)' 

READ( •, ' (Al ) ' )  ANS 

IF  (ANS  .EQ.  'N'  .OR.  ANS  . EQ .  'n')  GO  TO  30 
CALL  PPLOT(NSAnPLES,INFILE) 
C   •   Again? 

30   WRITE(*,*)'  ' 

WRITE  (LUNTI,*)  'Again?   (0-NO,  l-YES)  ' 
READ   (LUNTI,*)  lYESNO 
If  (  lYESNO.EQ.l  )  GO  TO  100 
C 

RETURN 

C   •   Foraats: 

700  FORKAT( 5X,F9.4,lX,f9.4,8X,I4) 
800  FORMAT( 5X,f9.4,lX,f9.4) 

90  0  FORMAT (//,2X, 'XPOS-' ,f6.2,2X, 'ZPOS-' , f 6 . 2 ) 

950  F0RMAT(2X,  'TEST  SECTION  VELOCITY- ', F5 . 1 , 2X , I  5 , '  SAMPLES') 
975  F0RMAT{2X, 'BLUE  AND  GREEN  MIXING  f REQUENCI ES-- ' , 2r5 . 1 , ' ( MHz ) ' ) 
END 
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SUBROUTINE  PT  PLOT( NSA«PLES , INFI LE 1 


*  Written  by  S.  D.  Hedrick  in  July,  1987  to  plot  computed  velocity  values 

*  corresponding  to  airfoil  phase  angles.   This  subroutine  is  intended  to 

*  provide  the  user  with  graphical  information  on  acquired  LDV  velocities 

*  and  assist  in  making  a  determination  regarding  the  adequacy  of  the  data 

*  obtained  at  the  current  probe  position. 

IMPLICIT  NONE 

INTEGER   NSAMPLES.I 

REAL*4    XVAL(4000) ,yVALU( 4000) ,yVALV( 4000) ,XC0NTR(4) ,yCONTR( 4) 

CHARACTER'IO  INTILE 

CHARACTER* 11  XLABEL 

CHARACTER* 15  YLABEL 

CHARACTER*19  TITLE 

CHARACTER* 1   DUMHY 

OPEN  (UNIT-ll,riLE-INriLE,STATUS-'OLD' ) 

00  10  I-1,NSAHPLES 

READ( 11,100)  YVALU( I ) ,yVALV( I ) ,XVAL( I ) 
XVAL( I )-XVAL( I )/10.0 
10   CONTINUE 
100   rORHAT( 5X , r9 . 4 , IX . r9 .  4 , 8X , r4 . 0 ) 

c 

C   Set  up  the  coordinate  axis 
C 

XLABEL- 'PHASE  ANGLE' 

YLABEL- 'VELOCITY  (ft/S)' 

TITLE- 'COMPUTED  VELOCITIES' 

XCONTR( 1 )-6.0 

XC0NTR(2 )-0.0  / 

XCONTR( 3)-360.0 

XCONTR(4)-45.0 

YCONTRd  )-6.0 

YCONTR( 2 )--200.0 

YCONTRI 3 )-600.0 

YCONTR( 4 )-100.0 

CALL  LGP$ PLOT ( 1, 'EXSY' ,XVAL,yVALU,0, XLABEL, YLABEL,  ,  ,  ,XCONTR, 
1  YCONTR,  , TITLE) 

C 

C      Plot    the    individual    data    points   on    the   axis 
C 

CALL    LGP$P0INT(1,XVAL,YVALU, .03, 3,NSAMPLES,     ,     ,     ) 

CALL    LGP$POINT( 1,XVAL,YVALV, .03,4,NSAMPLES,     ,     ,     ) 

CALL    LGP$PUT_TEXT(1,4.5,5.75, '•       U-COMPONENT' ,     ,     ) 

CALL    LGPSPUT_TEXT( 1,4. 5,5.5, '0      V-COHPONENT' ,     ,     ) 
C 

WFITE(*,*)    'Press    <CR>    to    terminate   plot' 

READ(*,200)    OUHMY 
200    PORMAT(Al) 

CALL    LGP$TERMINATE_PL0T(1,1) 

CLOSE(UNIT-ll) 

RETU1U4 

END 
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SUBROUTINE  PPLOT( NSArtPLES , INFI LE ) 
IMPLICIT  NONE 
INTEGER   NSAHPLES.I 

REAL       VMIN,V«AX,INTV 

REAL* 4    XVALI  4  000)  ,WALU(4000)  ,  YVALV(  4  000)  ,XCONTR(  4)  ,YC0NTR(4) 

CHARACTER'IO  INFILE 

CHARACTER*!!  XLABEL 

CHARACTER* 15  YLABEL 

CHARACTER*19  TITLE 

CHARACTER*!   DUMMY 

OPEN  (UNIT-!!,riLE-INFILE,STATUS-'OLD' ) 

DO  10  I-!,NSAMPLES 

READ( 11,100)  YVALU(l) ,YVALV(I) ,XVAL( I ) 
XVAL( I ) -XVALI I )/10.0 
10   CONTINUE 
100   FORnAT(5X,F9.4,lX,r9.4,8X,r4.0) 
C 

C   Set  up  the  coordinate  axis 
C 

WRITE( *,*) 'ENTER  VELOCITY  RANGE  DESIRED  (Bin, max)' 

READf*,*)  VMIN,VMAX 

INTV-(VMAX-VMIN)/5. 

XLABEL- 'PHASE  ANGLE' 

YLABEL- 'VELOCITY  (ft/B)' 

TITLE- 'COMPUTED  VELOCITIES' 

XCONTRI 1 )-6.0 

XCONTR(2)-0.0 

XCONTR( 3)-360.0 

XCONTRI 4 )-45.0 

YCONTR( 1 )-6.0 

YCONTRI 2 )-VMIN 

YCONTR( 3)-VMAX 

YCONTRI 4 )-INTV 

CALL  LGP$PLOT( 1 , 'EXSY' ,XVAL,YVALU,0, XLABEL, YLABEL,  ,  ,  ,XCONTR, 
1  YCONTR,  , TITLE) 

C 

C   Plot  the  individual  data  points  on  the  axis 
C 

CALL  LGP$POINT( 1,XVAL,YVALU, .03,3,NSAMPLES,  ,  ,  ) 

CALL  LGPSPOINTI l.XVAL, YVALV, .0  3,4,NSAMPLES,  ,  ,  ) 

CALL  LGP$PUT_TEXT(1, .05, . 3, 'ASTERISK-U  COMPONENT',  ,  ) 

CALL  LGPSPUT  TEXT( 1 ,. 05 ,. 1 ,' CIRCLE-V  .COMPONENT' ,  ,  ) 
C 

WRITE(*,*)  'Press  <CR>  to  terminate  plot' 

READ(*,200)  DUMMY 
200  FORMAT! Al) 

CALL  LGP$TERMINATE_PL0T(1,1) 

CLOSE(UNIT-ll) 

RETURN 

END 
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SUBROUTINE  HIST( YU , YV, YMAXU , YNAXV ) 


•  Written  by  S.  D.  Hedrick  in  August,  1987  to  plot  a  histogram  of  the  per- 

•  centage  of  data  points  which  fall  within  specified  velocity  limits.  The 

•  subroutine  produces  histograms  for  each  of  two  velocity  components.   Scal- 

•  ing  parameters  for  the  vertical  axis  and  data  arrays  are  computed  within 

•  a  preceding  subroutine  (HWWA). 
* 

IMPLICIT  NONE 

REAL '4  YU(160),yV(120) ,XC0NTRU(4) ,XC0NTRV(4) ,YC0NTRU(4) 

REAL '4  YUAXU , YHAXV , YCONTRV ( 4 ) 

REAL  XLOW,XHI,XLOV«ru(  160),XHIU(160)  .XL0WV(  120),XHIV(120) 

INTEGER  I 

CHARACTER*29  UTITLE 

CHARACTER*27  VTITLE 

CHAitACTER*10  XLABEL 

CHARACTER* 18  YLABEL 

CHARACTER* 1  DUHHY 
C 

C   Setting  up  the  coordinate  axis 
C 

XLABEL- ' VELOCITIES ' 

YLABEL- 'PERCENTAGE  OF  DATA' 

UTITLE- 'HORIZONTAL  VELOCITY  HISTOGRAM' 

VTITLE-'VERTICAL  VELOCITY  HISTOGRAM' 

XCONTRU( 1 )-6.0 

XCONTRU( 2 )--200. 0 

XCONTRU( 3)-600.0 

XCONTRU( 4 )-100.0 

XCONTRV( 1 )-6.0 

XCONTRVI 2)--300.0 

XCONTRV( 3)-300.0 

XCONTRV( 4 )-100.0 

YCONTRU( 1 )-6.0 

YCONTRU( 2) -0.0 

YCONTRU( 3 )-YMAXU 

YCONTRU( 4 ) -10.0 

YCONTRVI 1 )-6.0 

YCONTRV( 2 )-0.0 

YCONTRV( 3)-YMAXV 

YCONTRVI 4 ) -10.0 

CALL  LGP$PLOT( 1 , ' IXSY' ,  , YU , 0 , XLABEL , YLABEL ,  ,  ,  , XCONTRU , 
1  YCONTRU,  , UTITLE) 


c 

c 

Setting  up  arrays 

c 

XLOW--200.0 
XHI  — 195.0 
DO  10  1-1,160 
XLOWU(I)-XLOI» 
XHIU( I )-XHI 
XLOI»-XLOW*5.0 
Xai-XHI-t-5.0 

10 

CONTINUE 
XLOW— 300.0 
XHI  — 295.0 
DO  20  X-1,120 
XL(WV(I)-XLOM 

for  horizontal  and  vertical  histogram  bar  parameters 
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XHIV( I )-XHI 
XLOW-XLOWt-5.0 
XHI-XHI-t-5.  0 
20   CONTINUE 
C 

C   Plotting  the  horizontal  velocity  histogram 
C 

CALL  LGP$HIST(1,XLOWU,XHIU,VU,160,  ,  ,2) 
C 

WRITE! *,*)  'PRESS  <CR>  TO  TERMINATE  PLOT' 
READ! *, ' (Al) ' )  DUMMY 
CALL  LGPSTERMINATE_PL0T(1,1) 
C 

C   Repeating  the  procedures  to  plot  the  vertical  velocity  histograi 
C 

CALL  LGP$PLOT( 1, ' IXSY' ,  , W, 0 , XLABEL , YLABEL ,  ,  ,  ,XCONTRV, 
1  YCONTRV,  .VTITLE) 

CALL  LGP$HIST( 1 ,XLOWV,XHIV,YV,120,  ,  ,2) 
C 

WRITE(*,*)  'PRESS  <CR>  TO  TERMINATE  PLOT' 

READ! • , ' (Al ) ' )  DUMMY 

CALL    LGP$TERMINATE_PLOT( 1,1) 

RETURN 

END 
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PROGRAM  ORGIZE 


Written  by  S.  D.  Hedrick  in  July,  1987  to  sequentially  organize  data  files 
to  obtain  a  least  value  of  phase  angle  (third  column  of  data),  and  associ- 
ated velocities,  at  top  of  file.   This  rearrangement  of  data  was  desired 
to  assist  with  subsequent  manipulation  and  analysis.   In  September,  1987 
this  routine  was  incorporated  as  the  central  program  in  the  processing  and 
analysis  of  unsteady  LDV  velocity  data  generated  by  an  oscillating  airfoil. 
The  program  calls  assisting  subroutines  to  disgard  bad  data  points  and  pro- 
vide wake  profile  plots. 


IMPLICIT  NONE 

REAL  X( 4000 ),Y( 4000) , XPOS , ZPOS , VELO , UVL ,  \A:L , UNDIM( 30) ,VNDIM( 30) 

REAL  2P0( 30) , SAVl , SAV2 , MENU , MENV, PREJ , MEANU ( 18) ,MEANV( 18) 

REAL  MPHAVG( 18) 

INTEGER  Z( 4000) , NSAMPLES , PHASA, I , J , K , N , H, SAVE , REJ , STDY 

CHARACTER'IO  INTILE 

CHARACTER* 1  ANS 

STDY-0 

•)'IS  ANALYSIS  rOR  NON-OSCILLATING  AIRFOIL?  ( Y/N ) ' 
(Al ) ' )  ANS 

.EO-  'Y'  .OR.  ANS  .EQ.  'y')  STDY-1 

*) 'ENTER  NUMBER  OF  DATA  ACQUISITION  POINTS  IN  PROFILE', 
'  ( 30  maximum ) ' 
N 

'ENTER  PHASE  ANGLE  DESIRED  FOR  ANALYSIS  (degrees)' 
PHASA 


WRITE( • 
READ( *, 
IF  (ANS 
WRITE( • 


•) 


10 


READ( 

WRITE( * , • 

READ( *, • ) 

DO  70  K-1 

WRITE( • 

READ( », 

WRITE( • 

READ( * , 


,N 

,  •  )  'ENTER  NAME  OF 
' (AlO) ' )  INFILE 
,  •)  'ENTER  NUMBER 
'  )  NSAMPLES 


INPUT  DATA  FILE  ( LDVXXX . DAT ] 


OF  SAMPLES  IN  DATA  FILE' 


0PEN(UNIT-11 , FILE- INFILE, STATUS- 'OLD' ) 
DO  10  I-l, NSAMPLES 

READ(11,100)  X( I ) ,Y( I)  ,Z(I) 
CONTINUE 

READ(11,200)  XPOS.ZPOS 
READ! 11, 300)  VELO 
IF  (STDY  .EQ.  1)  THEN 
CLOSE(UNIT-H  ) 

CALL  STEADY( NSAMPLES, INFI LE , MENU , MENV ) 
UNDIM( K )-MENU/VELO 
VNDIMI K )-MENV/VELO 
GO  TO  80 
END  IF 

M-NSAMPLES-1 
DO  20  J-1, NSAMPLES 
DO  30  I-1,H 

IF  (Z(I+1)  .LT.  Z(I))  THEN 
SAVE-Z(I) 
SAVl-X(I) 
SAV2-Y(I) 
Z(  D-Zd-t-l) 
X(I)-X(I+1) 
Y(I)-Y(lfl) 
Z(I-fl)-SAVE 
X(I+1)-SAV1 
Y( I*1)-SAV2 
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END  ir 
30  CONTINUE 
20     CONTINUE 

CLOSE(UNIT-ll) 

OPEN ( UNI T- 1 5, FILE- I NFILE, STATUS- 'OLD' ) 

DO  40  I-1,NSAMPLES 

WRITE! 15,100)  X( I ) ,Y( I ) ,Z(I ) 
40     CONTINUE 

CL0SE(UNIT-15) 

WRITE( •, •) 'DO  YOU  DESIRE  A  VELOCITY/PHASE  PLOT  OF 
1  '  DATA?   (Y/N — DEF-Y)' 

READ( *, ' (Al ) ' )  ANS 

IF  (ANS  .EQ.  'N'  .OR.  ANS  .EQ.  'n')  GO  TO  50 
CALL  PT_PLOT(NSA«PLES, INFILE) 
50     CALL  BDPTS(NSA«PLES, INFILE, MENU, MENV,REJ,PREJ, MEA; 
1  MPHAVG) 

CALL  AVG(NSAMPLES, INFILE) 

WRITE! *,*) 'DO  YOU  DESIRE  INTERPOLATION  OF  THE  DAT; 
READ( *, ' (Al ) • )  ANS 

IF  (ANS  .EQ.  'N'  .OR.  ANS  .EQ.  'n')  GO  TO  55 
CALL  INTERP(NSAnPLES, INFILE) 
55     CALL  CVFIT(NSA«PLES, INFILE, PHASA,UVL,WL, MENU, MENN 
1  KEANU , NEANV , MPHAVG ) 

UNDIM( K )-UVL/VELO 
VNDIM(  K)-WL/VELO 
80     ZPO(K)-ZPOS 
70   CONTINUE 

CALL  PR0PL0T(UNDIM,VNDIM,2P0,XP0S,N) 

:   Format  statements 

100   TORMAT( 5X,F9. 4,1X,F9.4,8X, 14) 
200   F0RMAT(//,7X, F6.2,7X,r6.2) 
300   FORMAT( 24X,F5.1 ) 

STOP 

END 
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SUBROUTINE  BDPTS { NSArtPLES , INTI LE , MENU , MENV , RE J , PRE J , 
1  «EANU,MEAhA:,MPHAVG) 

REAL  X(4000),Y(4000) ,«EANU( 18 ) ,MEANV( 18 ) ,STDVU( 18) 
REAL  STDW(  18  )  ,MPHAVG(  18)  ,  MINU  ,  MINV  ,  MAXU  ,  MAXV  ,  MENU  ,  MENV 
INTEGER  Z( 4000) , PHASEA, PHASES , SUM , COUNT , NSAMPLES , NSAMP , REJ 
REAL»4  XCONTRf 4 ) ,YCONTR( 4 ) ,PREJ 
CHARACTER*!  ANS , DUMMY , DUMMI 
CHARACTER'IO  INFILE 
CHARACTER*11  XLABEL 
CHARACTER*21  YLABEU 
CHARACTER*19  YLABEV 
CHARACTER*18  TITLE 

OPEN (UNIT- 15, FILE- INFILE, STATUS- 'OLD' ) 
NSAMP-NSANPLES 
DO  10  I-1,NSAMPLES 

READ(15,100)  X(  I  )  ,y( I ) ,Z( I ) 
10   CONTINUE 
100   rORMATC 5X,r9.4, lX,r9.4,8X,I4) 
PHASEA-0 
PHASEB-200 
DO  20  1-1,18 
SUM-0 
SUMU-0. 
SUMV-0. 
PHAVG-0. 
DO  30  J-1,NSAMPLES 

IF  (Z(J)  .LE.  PHASES  .AND.  2(J)  . GE .  PHASEA)  THEN 
SUMU-SUMUt-X(  J) 
SUMV-SUMV+Y( J) 

PHAVG-PHAVG+( FL0AT(2( J) )/10. ) 
SUM-SUM>1 
END  IF 
30     CONTINUE 

IF  (SUM  .EQ.  0)  THEN 
MEANUI I )-0. 
MEANV( I)-0. 

MPHAVG( I )-( FLOAT( PHASEA)* 100. )/10. 
VARU-0. 
VARV-0. 
GO  TO  4  0 
END  IF 

MEANU( I )-SUMU/FLOAT(SUM) 
MEANV( I )-SUMV/FLOAT(SUM) 
MPHAVG( I )-PHAVG/FLOAT(SUM) 
VARU-0. 
VARV-0. 
DO  40  J-1,NSAMPLES 

IF  (Z(J)  .LE.  PHASES  .AND.  Z(J)  . GE .  PHASEA)  THEN 
Al-X( J)-HEANU( I ) 
A2-Y( J)-MEANV( 1 ) 
B1-A1*«2/FL0AT(SUM) 
B2-A2**2/FLOAT(SUM) 
VARU-VARU-fBl 
VARV-VARV+B2 
END  IF 
40     CONTINUE 

STDVU( I )-SORT( VARU) 
STDWI  I  )-SQRT(VARV) 
PHASEA-PHASEAf200 
PHASBB-PBASEB4-200 
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20   CONTINUE 

WRITE( », * ) 'DO  YOU  DESIRE  A  STANDARD  DEVIATION  PLOT?  ( Y/N  -  DEF-Y) 
READ(*,200)  ANS 
2  00   rORMAT(Al) 

IF  (ANS  .EQ.  'N'  .OR.  ANS  . EQ .  'n')  THEN 

GO  TO  50 
END  IF 
C 

C   Setting  up  the  coordinate  axis  for  the  plot 
C 

XLABEL-' PHASE  ANGLE' 
YLABEU-' HORIZONTAL  VELOCITIES' 
YLABEV-' VERTICAL  VELOCITIES' 
TITLE- 'STANDARD  DEVIATION' 
XCONTR( 1 )-6.0 
XCONTR(2)-0.0 
XCONTR( 3)-360.0 
XCONTRI 4) -4 5.0 
YCONTRd  )-6.0 
yCONTR( 2)--200.0 
YCONTR( 3)-600.0 
YCONTRI 4 ) -100.0 

CALL  LGP$PLOT( 1 , 'EXSY' ,«PHAVG,HEANU,0, XLABEL, YLABEU,  ,  ,  , XCONTR , 
1  YCONTR,  .TITLE) 

C 

C   Plotting  the  horizontal  velocity  data  points  on  the  axis 
C 

CALL  LGP$POINT( 1,MPHAVG,HEANU, .03,4,18,  ,  ,  ) 
C 

C   Plotting  the  horizontal  standard  deviation 
C 

CALL  LGP$STNDEV( l,nPHAVG,HEANU,STDVU,18,  ) 
C 

45   WRITE(*,*)  'PRESS  <CR>  TO  TERMINATE  PLOT' 
READ(*,200)  DUMMY 
CALL  LGPSTERMINATE_PLOT( 1,1 ) 
C 

C   Plotting  the  vertical  standard  deviation 
C 

CALL  LGP$PLOT( 1, 'EXSY' ,«PHAVG,MEANV,0, XLABEL, YLABEV,  ,  ,  , XCONTR, 
1  YCONTR,  , TITLE) 

CALL  LGP$POINT( 1,MPHAVG,MEANV,  .03,  4,18,  ,  ,  ) 
CALL  LGP$STNDEV(  1  ,MPHAVG,MEANV,STDW,18,   ) 
47   WRITE(*,*)  'PRESS  <CR>  TO  TERMINATE  PLOT' 
READ(*,200)  DUMMI 
CALL  LGP$TERMINATE_PL0T(1,1) 
C 

50   COUNT- 1 
PHASEA-0 
PHASEB-200 
DO  60  1-1,18 

MINU-MEANUd  )-2.0*STDVU(  I  ) 
«INV-MEANV(  I  )-2.0*STDW(  I  ) 
HAXU-MEANUd  ) +2  .  0*STDVU(  I  ) 
HAXV-«EANV(  I  )*2.0*STDW{  I  ) 
DO  70  J-1,NSAMPLES 

ir  (Z(J)  .LE.  PHASEB  .AND.  Z(J)  .GE.  PRASEA)  THEN 
ir  (X(J)  .LT.  HINU  .OR.  X(J)  .GT.  MAXU)  THEN 
COUNT-COUNT* 1 
If  ( J  .EQ.  1)  GO  TO  70 
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DO  80  K-0,J-2 

X( J-K)-X( J-K-1) 
Y( J-K)-Y( J-F-1 ) 
2( J-K)-2( J-K-1 ) 
80  CONTINUE 

GO  TO  7  0 
END  ir 

IF  (Y(J)  .LT.  MINV  .OR.  Y(J)   . GT .  MAXV )  THEN 
COUNT-COUNT-f  1 
ir  (J  .EQ.  1 )  GO  TO  7  0 
DO  85  K-0,J-2 

X(J-K)-X( J-K-1 ) 
Y(J-K)-Y( J-K-1) 
Z(J-K)-Z( J-K-1 ) 
85  CONTINUE 

END  IF 
END  IF 
70     CONTINUE 

PHASEA-PHASEA+2  00 
PHASEB-PHASEB-^2  00 
60   CONTINUE 

CL0SE(UNIT-15) 

0PEN(UNIT-12, FILE- I NFILE, STATUS- 'OLD' ) 

DO  90  I-COUNT,NSAMPLES 

WRITE(12,100)  X(I ) ,Y( I)  ,Z(I) 
9  0   CONTINUE 

NSAMPLES-NSAMPLES+1 -COUNT 

CL0SE(UNIT-12 ) 

OPEN ( UNI T- 1 3, FILE-INFILE, STATUS- 'OLD' ) 

SUU-0. 

SUV-0. 

DO  95  I-1,NSAMPLES 

READ( 13,100 )  X( I ) , Y( I ) ,Z( I  ) 
SUU-SUU+X( I ) 
SUV-SUV+Y( I ) 
95   CONTINUE 

MENU-SUU/FLOAT( NSAMPLES ) 
MENV-SUV/FLOAT( NSAMPLES ) 
CL0SE(UNIT-13 ) 
REJ-NSAMP-NSAMPLES 

PREJ-( l.-( FLOAT (NSAMPLES) /FLOAT (NSAMP) ) )*100. 
PHASEA-0 
PHASEB-200 
DO  96  1-1,18 
SUM-0 
SUHU-0. 
SUMV-0. 
PHAVG-0. 
DO  97  J-1, NSAMPLES 

IF  (Z(J)   .LE.  PHASES  .AND.  Z(J)   . GE .  PHASEA)  THEN 
SUMU-SUMU*X( J) 
SUHV-SUMV+Y( J) 

PHAVG- PHAVG+ ( FLOAT ( Z ( J )  ) /1 0 .  ) 
SUM-SUH-fl 
END  IF 

97  CONTINUE 

IF  (SUH  .EQ.  0)  THEN 

HeANU( I )-0. 

HEANVI I )-0. 

nPliAVG(  I  )-(FLOAT(  PHASEA) +  100.  )/10. 

GO  TO  98 
END  IP 

HEANU ( I ) -SUMU/FLOAT (  SUM ) 
MEANV( I )-SUMV/FLOAT(SUM) 
HPHAVG( I )-PHAVG/FLOAT(SUM) 

98  PHASEA-PHASEA'*'200 
PHASEB-PHASEB•^200 

96   CONTINUE 
RETURN 
END 
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SUBROUTINE  AVG( NSArtPLES , INFI LE ) 

REAL   X( 4000  )  ,Y( 4000) 

INTEGER   Z( 4000) , COUNT , SUM , NSAMPLES 

CHARACTER*10  INFILE 

COUNT-1 

M-NSAMPLES-1 

OPEN (UN IT- 12, FILE- INFILE, STATUS- 'OLD' ) 

DO    10    I-1,NSAI1PLES 

READ(12,100)  X( I) ,Y( I) ,Z( I) 
10   CONTINUE 
100   FORMAT( 5X,F9.4,lX,r9.4,8X,I4) 
DO  20  I - 1  ,  H 

IF    (2(1  +  1)     .EQ.     Z(I)     .AND.     Z(I>2)     .  NE .     Z(I))    THEN 
COUNT-COUNT* 1 
X(  1-^1)-(X(  I  +  D+Xd)  )/2.0 

y(i+i)-(y( i+i )*Y(i) )/2.0 

DO    30    J-0,I-2 

X( I-J)-X( I-J-1) 

Y(i-j)-y(i-j-i) 

Z( I-J)-Z(I-J-l) 
30       CONTINUE 
END  IF 

IF  (Z(I+2)  .EQ.  Z(I))  THEN 
SUM-1 
DO  25  J-2,NSAMPLES-I 

IF  (Z(I-fJ)  .EQ.  Z(I))  THEN 

SUM-SUMt-1 
END  IF 
25       CONTINUE 

DO  27  J-1,SUH 

X( I*SUM)-X( I*SUM)+X( I+SUM-J) 
Y(  It-SUH)-Y(  I+SUM)+Y{  I  +  SUM-J) 

27  CONTINUE 

X(  IfSUM)-X(  I-»-SUM)/(SUM+l) 
Y(  I+SUM)-Y(  IfSUM)/(SU«-t-l  ) 
DO  28  J-1,I-1 

X( I*SUM-J)-X( I-J) 

Y(  H-SUM-J)-y( I-J) 

2( I+SUM-J)-Z(I-J) 

28  CONTINUE 
COUNT-COUNT-^SUM 

END  IF 
20   CONTINUE 

CL0SE(UNIT-12) 

OPEN (UNIT- 14, FILE- INFILE, STATUS- 'OLD' ) 

DO  40  I-COUNT,NSA«PLES 

WRITE(14,100)  X( I) ,Y(I) ,Z( I ) 
4  0   CONTINUE 

NSAHPLES-NSAHPLES+l -COUNT 

CL0SE(UNIT-14) 

RETURN 

END 


64 


SUBROUTINE  INTERP ( NSArtPLES , INFILE) 

IMPLICIT  NONE 

CHARACTER*10   INPILE 

REAL   X( 3600) ,Y( 3600) 

INTEGER   Z( 36  00) .NSAMPLES, I , J , COUNT, NSA,  PAC 

CHARACTER*!   ANS 

OPEN (UNIT- 11 ,riLE-INFILE, STATUS- 'OLD' ) 

DO  10  I-l,NSArtPLES 

READ( 11,100)  X( I ) , Y( I ) ,Z( I) 
10   CONTINUE 
100   FORMAT( 5X,r9.4,lX,r9. 4,8X,I4) 
COUNT-0 
NSA-NSAMPLES 
DO  20  1-1,3600 

IF  (COUNT  .GE.  NSA)  THEN 
2( I  )-I-l 

X(I)-((X(1)-X(I-1))/(Z(1)-2(I-1)))*(Z(I)-Z(I-1))+X(I-1) 
y(I)-((Y(l)-Y(I-l))/(Z(l)-Z(I-l)))*(Z{I)-Z(I-l))+Y(I-l) 
NSAHPLES-NSAHPLES-*-! 
GO  TO  20 
END  IF 

IF  (Z(I)  .EQ.  I-l)  COUNT-COUNT+1 
IF  (2(1)  .NE.  I-l)  THEN 
NSArtPLES-NSAMPLES+1 
DO  30  J-0,NSA«PLES-I-1 

X(NSAJ1PLES-J  )-X(NSAMPLES-J-l) 
Y(NSAMPLES-J )-Y(NSA«PLES-J-l) 
Z(NSA«PLES-J)-2(NSArtPLES-J-l ) 
30        CONTINUE 
2(  I  )-I-l 
IF  (I  .EQ.  1)  THEN 

X(I)-((X(I*1 )-X(NSAMPLES) )/(2( I+l ) -Z ( NSArtPLES ) ) )*(Z( I )- 
1  Z  (  NSAMPLES  )  )  ■t-X(  NSAMPLES  ) 

Y{I)-((y(I*l )-Y(NSArtPLES) )/(Z( I+l) -Z ( NSAMPLES ) ) )*(2( I)- 
1  2  (  NSAMPLES  )  )-fY(  NSAMPLES) 

GO  TO  20 
END  IF 

X(I)-((X(I  +  l)-X(I-l))/(2(I  +  l)-Z(I-l)))*(Z(I)-Z(I-l))-t-X(I-l) 
Y(I)-((Y(1  +  1)-Y(I-1))/(2(I*1)-2(I-1)))M2(I)-2(I-1))*Y(I-1) 
END  IF 
20   CONTINUE 

CL0SE(UNIT-11 ) 
C       WRITE( *, • ) 'DO  YOU  DESIRE  TO  REDUCE  THE  NUMBER  OF  INTERPOLATED' 
C      1  '  DATA  POINTS  BY  AN  INTEGER  FACTOR  OF  36007  ( Y/N ) ' 

C       READ( *, ' (Al ) ' )  ANS 

OPEN (UNIT- 12, FILE- INFILE, STATUS- 'OLD' ) 
C       IF  (ANS  .EQ.  'Y'  .OR.  ANS  . EQ .  'y')  THEN 
C  WRITE (»,») 'ENTER  INTEGER  FACTOR' 

C  READ(*,*)  FAC 

FAC-20 

DO  3  4  I -FAC, NSAMPLES, FAC 

WRITE( 12,100)  X(I ) ,Y( I ) ,Z(I) 
34     CONTINUE 

NSAMPLES- 36 OO/FAC 
C  GO  TO  40 

C       END  IP 
C       DO  40  I-I, NSAMPLES 
C         WRITE(12,100)  X(I),Y(I),Z(I) 
C   4  0   CONTINUE 

CL0Se(UNIT-12) 

RETURN 

END 
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SUBROUTINE  FOUHl ER( F , M , A, B , N ) 

DIMENSION  F( 0:«), FA ( 0: 300 ),rB( 0:300) 

DATA  PI/3.14159265/ 

H-2*PI/« 

TYPE*,  H 

J-0 

DO  2  I-0,M 

C-4 

IF  (I  .EQ.  J)  C-2 

IF  (I  .EQ.  J)  J-J-^2 

FA( I )-F( I )*COS(N*H*I ) 

FBI  I )-F( I )*SIN(N*H*I ) 

IF  (I  .EQ.  M)  T3-FA(I) 

IF  (I  .EQ.  M)  TT3-FB(I) 

IF  (I  .EQ.  M-1)  T2-FA(I) 

IF  (I  .EQ.  «-l)  TT2-FB(I) 

IF  (I  .EQ.  M-2)  Tl-FA(I) 

IF  (I  .EQ.  M-2)  TTl-FB(I) 

IF  (I  .EQ.  M-3)  TO-FA(I) 

IF  (I  .EQ.  M-3)  TTO-FB(I) 

I F  (  I  . EQ .  0)  GO  TO  2 

FA(  I  )-FA(  I-l  )-t-C*FA(I) 

FB( I )-FB( I-1)*C*FB(1) 
CONTINUE 

A-H/3./PI*( rA(M)-T3) 
B-H/3./PI*( FB(M)-TT3) 

IF  (J  .NE.  M-t-1  )  GO  TO  3 

A-H/3./PI*( FA(M)-4. •T3-2 . *T2-4 . •Tl-TO) 

A-A+3.  *H/8./PI*(T0+3.  ♦Tl-t-3.*T2  +  T3)- 

B-H/3./PI*( FB(M)-4 .•TT3-2. •TT2-4. 'TTl-TTO) 

B-Bt-3.  •H/8./PI*(TT0+3.*TTl  +  3.  •TT2  +  TT3) 

CONTINUE 

RETURN 

END 
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SUBROUTINE    CVriT(  NSArtPLES  ,  INFILE,  PHASA  ,  UVL  ,  WL  ,  MENU  ,  MENTV,  REJ  ,  PREJ  , 
1  «EANU,MEANV,MPHAVG) 

» 

»  Written  by  S.  D.  Hedrick  in  October,  1987  to  fit  a  sinusoidal  curve  to  LDV 

'  velocity  data  generated  by  an  oscillating  airfoil.   The  subroutine,  utiliz- 

'  ing  FOURIER,  adapted  from  June  1977  OCS  Airfoil  program,  fits  a  single 

*  cycle  sine  wave  to  the  data  and  computes  the  u-  and  v-velocity  components 

*  at  a  phase  angle  specified  for  the  analysis. 
» 

IMPLICIT  NONE 
CHARACTER* 1  ANS , DUMMY 
CHARACTER*10  INFILE 

INTEGER  NSAflPLES,  I  ,  PHASA  ,  NSA  ,  REJ  ,  ARG 

REAL  FU( 0 : 3600 ), rv ( 0: 3600 ),AU,AV,BU,BV, PI, X( 3600 ),Y( 3600) 
REAL  UVL,WL,MENU,MENV,PREJ,MEANU(  18)  ,MEANV(  18)  ,MPHAVG{  18) 
DATA  PI/3.1415927/ 

WRITE( *,•) 'DESIRE  A  VELOCITY/PHASE  PLOT  OF  THE  REFINED  DATA?' 
WRITE! »,*)'( Y/N — DEF-Y) ' 
READ! • , ' (Al ) ' )  ANS 

IF  (ANS  .EQ.  'N'  .OR.  ANS  . EQ .  'n')  GO  TO  7 
CALL  RPLOT( NSAMPLES , INFILE, RE J, PRE J ) 
7   OPEN (UN IT- 11 , FILE- INFILE, STATUS- 'OLD' ) 
DO  5  I-l ,NSA«PLES 

READ( 11,100)  X(I )  ,Y(I  ) 
5   CONTINUE 
100   FORMAT( 5X, F9.4,1X,F9.4) 
CL0SE(UNIT-11  ) 
DO  10  I-l, NSAMPLES 
FU(  I-l )-X(  I  ) 
FV(  I-l )-Y( I  ) 
10   CONTINUE 

NSA-NSAMPLES-1 

CALL  FOURIER( FU,NSA,AU,BU,1) 
CALL  FOURIER! FV, NSA, AV,BV,1 ) 
DO  20  I-0,NSA 

FU(  I  +  l  )-MENU*AU*COS(  (  2.  •  PI/FLOAT (  NSAMPLES  )  )* FLOAT (  I  )  )  ••• 
1  BU»SIN( ( 2. "PI/FLOAT! NSAMPLES) )• FLOAT ( I ) ) 

FV(  I  +  l  )-MENV*AV*COS(  (  2.  'PI/FLOAT!  NSAMPLES  )  ) 'FLOAT!  I  )  )  ••• 
1  BV*SIN( ( 2. 'PI/FLOAT (NSAMPLES) )• FLOAT ( I ) ) 

20   CONTINUE 

WRITE( •,*)' PLOTS  WILL  NOW  BE  PROVIDED  TO  COMPARE  SINUSOID  WITH', 
1  '  VELOCITY  DATA.    PRESS  <CR>  TO  CONTINUE.' 

READ! • , ' (Al ) ' )  DUMMY 

CALL  PLOT ( NSAMPLES , INFILE , fU , FV , MENU , MENV , MEANU , HEANV, 
1  MPHAVG) 

ARG-NSAMPLES*(FLOAT(PHASA)/360. ) 
WRITE( • , * ) '  ARG  -' ,ARG 
UVL-FU(ARG) 
WL-FV!ARG) 
RETURN 
END 
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SUBROUTINE  PLOT( NSAHPLES , INF I LE , FU , TV , MENU , MENV , MEANU , 
1  HEANV.HPHAVG) 

» 

»   Written  by  S.  D.  Hedrick  in  October,  1987  to  plot  a  comparison  of  velocity 

*      data  and  a  sinusoidal  curve  fit. 

t 

IMPLICIT  NONE 

INTEGER   NSAMPLES, I , BOUND , STOPU , STOPV 

REAL '4    XVAL( 4  000) , YVALU( 4  0  00) ,yVALV( 4000) ,XCONTR( 4 ) , YCONTU( 4) 

REAL •4    yCONTV( 4 ) , FU(NSAMPLES) , FV(NSAMPLES) 

REAL • 4    YM I NU , YMAXU , YM I NV , YMAXV , M  ENU , MENV 

REAL       MEANU ( 18) ,MEANV( 18 ) ,MPHAVG( 18) 

CHARACTER*10  INFILE 

CHARACTER*11  XLABEL 

CHARACTER*15  YLABEL 

CHARACTER*21  UTITLE 

CHARACTER*19  VTITLE 

CHARACTER* 1   DUMMY 

OPEN (UNIT- 11, FILE- INFILE, STATUS- 'OLD' ) 

DO  10  I-1,NSAMPLES 

READ( 11 , 100)  YVALU( I ) ,YVALV( I ) , XVAL ( I ) 
XVAL( I )-XVAL( I )/10. 
10   CONTINUE 
100   FORMAT*  5X,F9.4, 1X,F9.4,8X,F4.0) 

XLABEL- 'PHASE  ANGLE' 

YLABEL- 'VELOCITY  (ft/S)' 

UTITLE- 'HORIZONTAL  VELOCITIES' 
VTITLE-'VERTICAL  VELOCITIES' 

YMINU-MENU-150. 

YMAXU-MENUt-150. 

YMINV-MENV-150. 
YMAXV-MENV4  150. 
XCONTR( l)-6.0 
XCONTR( 2) -0.0 
XCONTR( 3 )-360.0 
XCONTRI 4 )-45.0 
YCONTU( 1 )-6.0 
YCONTUI 2 )-YMINU 
YCONTUI 3)-YMAXU 
YCONTU( 4 )-50.0 
YCONTV( 1 )-6 .0 
YC0NTV(2 )-YMINV 
YCONTVI 3 )-YMAXV 
YCONTVI 4 )-50.0 

CALL  LGPSPLOTd, 'EXSY' , XVAL,  FU,NSAMPLES, XLABEL,  YLABEL,  ,  ,  , 
1  XCONTR,YCONTU,  , UTITLE) 

CALL  LGP$POINT( 1,MPHAVG, MEANU, .03, 3,18,  ,  ,  ) 
WRITE! •,*) 'PRESS  <CR>  TO  TERMINATE  PLOT' 
READ(*,2  00)  DUMMY 
CALL  LGP$TERMINATE_PLOT( 1 ,1 ) 

CALL  LGPSPLOTd, 'EXSY'  , XVAL, fV,NSAHPLES, XLABEL, YLABEL,  ,  ,  , 
1  XCONTR,YCONTV,  , VTITLE) 

CALL  LGP$P0INT(1,MPHAVG,MEANV, .03, 3, 18,  ,  ,  ) 
WRITB(»,*)  'Press  <CR>  to  terminate  plot' 
REAO(*,200)  DUMMY 
200  FORMAT(Al) 

CALL  LGP$TERMINATE_PL0T(1,1) 

CLOSe(UNIT-ll) 

RETURN 

END 
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SUBROUTINE  PROPLOT ( UNDIM , VNDIN , ZPO, XPOS , N ) 
IHPLICIT  NONE 

REAL* 4  UNDIM( 30) , VNDIW( 30 ) ,2P0( 30) ,XCONTU( 4 ) ,XCONTV( 4 ) 

REAL •4  YCONTR( 4 ) ,UNDI( 3000 ) ,VNDI( 3000 ),ZP( 3000) 

REAL  XPOS 

INTEGER    N,NBR 

CHARACTER»29    XLABEU , XLABEV 

CHARACTER*32  YLABEL 

CHARACTER»27  UTITLE 

CHARACTER*25  VTITLE 

CHARACTER*!   DUMMY, ANS 

CHARACTER*9       XPO 

CHARACTER*4   XXPO 

XLABEU-'VELOCITY  RATIO  ( UAJ  ( INFIN )  )  ' 

XLABEV- 'VELOCITY  RATIO  ( V/U( INFIN )) ' 

YLABEL-' PROBE  POSITION  (IN.  FROM  CENTER)' 

UTITLE-'HORIZONTAL  VELOCITY  PROFILE' 

VTITLE-'VERTICAL  VELOCITY  PROFILE' 

XPOS-XPOS/3.0 

WRITE(UNIT-XXPO,FMT-' ( F4 . 1 ) ' ) XPOS 

XPO-'X/C  -'//XXPO 

XCONTUI 1 )-6.0 

XCONTU( 2 )-0.0 

XCONTU( 3)-l .2 

XCONTU( 4 )-0.2 

XCONTV( 1 )-6.0 

XCONTV( 2)--0.7 

XCONTV( 3) -0.7 

XCONTVI 4 )-0 . 2   ' 

YCONTR( 1 )-6.0 

YCONTR( 2)--2.0 

YCONTR( 3 )-2.0 

YCONTR( 4 )-0.5 

NBR-10*(N-1 )*1 

CALL  LGP$PLOT( 1, 'EXSY' , UNDIM , ZPO , 0 , XLABEU , YLABEL , 
1  ,  ,  ,XCONTU, YCONTR,  , UTITLE) 

CALL  LGPSPOINT( 1 ,UNDIM,ZPO, .03, 3,N,  ,  ,  ) 

CALL  LGPSPUT_TEXT( 1 , .4, 5. 5,XP0,  ,  ) 

WRITE( *,•)' PRESS  <CR>  TO  TERMINATE  PLOT' 

READ{ *, ' (Al ) ' )  DUMMY 

CALL  LGPSTERMINATE_PLOT( 1 ,1 ) 

CALL  LGPSPLOT( 1, 'EXSY' , VNDIM , ZPO , 0 , XLABEV , YLABEL , 
1  ,  ,  ,XCONTV, YCONTR,  , VTITLE) 

CALL  LGPSPOINT{ 1 , VNDIM, ZPO, .03,3, N,  ,  ,  ) 

CALL  LGPSPUT_TEXT( 1,  .4,5.5,XPO,  ,  ) 

WRITE( •,»)' PRESS  <CR>  TO  TERMINATE  PLOT' 

READ{ *, ' (Al ) ' )  DUMMY 

CALL  LGPSTERMINATE_PLOT( 1, 1 ) 

WRITE( *, • ) 'DO  YOU  DESIRE  A  SPLINE  INTERPOLATION  CURVE  FIT  OF' 
1  '  THE  DATA  POINTS?  ( Y/N ) ' 

READ( •, ' (Al) ' )  ANS 

IF  (ANS  .EQ.  'N'  .OR.  ANS  . EQ .  'n')  GO  TO  10 

CALL  LGP$PLOT( 1 , 'EXSY' , UNDIM , ZPO , 0 , XLABEU , YLABEL, 
1  ,  ,  ,XCONTU, YCONTR,  .UTITLE) 

CALL  LGP$POINT( 1,UNDIM,ZP0, .03, 3,N,  ,  ,  ) 

CALL  LGPSSPLINE  QHC ( UNDIM , ZPO , N , 1 0 , UNDI , ZP , ) 

CALL  LGPSPLOTMd,  'EXSY'  ,UNDI  ,ZP,NBR,  ) 

CALL  LGPSPUT_TEXT(1, .4,5.5,XPO,  ,  ) 
WRITE( •,*) 'PRESS  <CR>  TO  TERMINATE  PLOT' 

READ( *, ' (Al ) ' )  DUMMY 
CALL  LGP$TERMINATE_PLOT( 1, 1 ) 

CALL  LGP$PLOT( 1 , 'EXSY' , VNDIM , ZPO , 0 , XLABEV, YLABEL, 
1  ,  ,  ,XCONTV, YCONTR,  , VTITLE) 

CALL  LGPSPOINT( 1,VNDIM,ZP0. .03, 3,N,  ,  ,  ) 
CALL  LGPSSPLINE  QHC ( VNDIM , ZPO , N , 1 0 , VNDI , ZP , ) 
CALL  LGPSPLOTMd,  'EXSY'  ,VNDI  ,ZP,NBR,  ) 
CALL  LGP$PUT_TEXT(1 , .4,5.5,XPO,  ,  ) 
WRITE( •,*) 'PRESS  <CR>  TO  TERMINATE  PLOT' 
READ( ♦, ' (Al ) ' )  DUMMY 
CALL  LGPSTERMINATE  PL0T(1,1) 
10   RETURN  ~ 

END 
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